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Py 

Melville Dannehl Aldridge, III^, Gary E. Copeland^, 
and Charles N. Harward^ 

ABSTRACT 

Measurements of gaseous atmospheric ammonia from 
September 1980 to mid-March 1981 were obtained with an 
experimental system employing laser photo-acoustic 
spectroscopic techniques. Ammonia concentrations reached 
minimal levels (approximately Q.I ppb) in early winter, 
followed by a sudden late winter increase. A direct 
relationship between ambient ammonia levels and air 

temperature was inferred from the data (linear correlation 
coefficient r«0.53). Ammonia concentrations were 
determined to be directly related to the absolute humidity 
of the air (r«0.72}; a weaker relationship between ammonia 
concentrations and relative humidity was discovered 
(r«0.37). The data also indicated that ammonia levels were 
generally higher within continental air masses than those 
of maritime origin. Soil parameters such as pH and 

moisture content were found to have a major bearing on the 
release of gaseous ammonia from soils in the region. 

^Graduate Research Assistant, ^Research Assistant Professor, 
and ^Associate Professor, Department of Physics, Old Dominion 
University, No’-folk, Virginia 23508. 
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INTRODUCTION 


Purpose 

In July 1980 a grant from the National Aeronautics and 
Space Administration (NASA) was awarded to Old Dominion 
University (principal investigator: Dr. G.B. Copeland) 
to conduct research into the temporal distribution of 
gaseous ammonia in the air of Southeastern Virginia, Since 
the late summer of 1980 measuremente of ambient ammonia 
have been taken with an experimental measurement system 
housed at NASA Langley Research Center in Hampton, 
Virginia. This paper will present those measurements and 
will investigate various factors which may have influenced 
the variations in ammonia concentrations noted during the 
course of this research. An attempt will also be made to 
place the measurements in a proper perspective. That is, 
results of a thorough review of the current knowledge 
concerning atmospheric ammonia will be presented. The 
review will contain information on the role of ammonia in 
the nitrogen cycle, aspects of ammonia’s atmospheric 
chemistry that may have deleterious consequences in the 
environment, and the results of other research programs 
which included measurements of atmospheric ammonia. In 
addition, the measurement technique used for this study 
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will be described in some detail— ‘and other methods of 
monitoring atmospheric ammonia will be discussed-- 0 o that 
the existing technologies can compared with the 
experimental technique employed in this study. It has been 
recognized for some time that gaseous ammonia is an 
important trace constituent of the atinosphere ; yet accurate 
measurement of it has proven to be extraordinarily 
difficult. As a result, relatively little is now known 
about ammonia in the atmosphere. The primary purpose of 
this peper will be to increase our knowledge and help pave 
the way for future research in this area. 

Atmospheric Ammonii* and the Nitrogen Cycle 

Gaseous ammonia, the familiar compound consisting ox; a 
nitrogen atom bound to three hydrogen atomt , is a 
ubiquitous trace constituent of the atmosphere. Although 
poisonous in high concentrations, gaseous ammonia rarely 
constitutes a health hazard outside of the laboratory. In 
the environment the ammonia molecule, together with the 
ammonium ion (NH4+), represent the most common forms of 
fixed nitrogen, which is essential to life {!]. Ammonia is 
therefore generally considered a beneficial substance. 
Ammonia is so essential for food production that in the 
last fifty years man has increasingly supplemented the 
enormous amount of ammonia produced by natural processes 
with manufactured ammonia. He has also planted large 
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numbers of nitrogen-fixing plants to augment natural 
production of fixed nitrogen. Man's activities account for 
a large and increasing share of the totalj quantity of fixed 
nitrogen produced each year. Although estimates vary 
widely# by most accounts the quantity of nitrogen fixed 
annually by industry and legume crops approximately equals 
that fixed "naturally" (1). How much of this fixed 
nitrogen actually finds its way into the atmosphere as 
gaseous ammonia is not known. Better understanding of the 
nitrogen cycle, and especially of soil processes, is 
necessary before scientists can assay what effect the 
increase in fixed nitrogen production will have on the 
atmospheric pool of ammonia. 

Natural production of ammonia occurs as an integral 
part of the nitrogen cycle. Organic nitrogen in soil and 
in water is mineralized by microorganisms into ammonium or 
nitrate ions, forms of nitrogen suitable for assimilation 
by Alants. In plant tissues the nitrogen is converted to 
organic form again, only to be returned to the soil or 
water by death and decay. However, the cycle. is 
complicated by the action of some species of microbes, 
which extract energy by converting nitrate to gaseous 
nitrogen or nitrous oxide. This process, called 
denitrification, must be balanced by an input of fixed 
nitrogen for the cycle to remain in equilibrium. Most 
fixed nitrogen in the soil is provided, by certain species 
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of microorganisms, usually existing in symbiotic 
relationships with plants, which metajol ically convert 
atmospheric nitrogen to ammonium. A smaller amount o£ 
fixed nitrogen is created in the atmosphere by ionizing 
processes such as lightning. 

The processes which lead to denitrification are a 
major sink for fixed nitrogen in the soil. Moreover, 
volatilization of ammonia from the soil into the atmosphere 
also constitutes a large sink for fixed nitrogen in the 
soil and represents the major source of atmospheric 
ammonia. According to Dawson [8) some 50 million metric 
tons per year of ammonia are released by volatilization 
from the soil. Once in the atmosphere gaseous ammonia 
resides there for comparably short periods. Research 
indicates that the average residence time is approximately 
five to ten days [1] [271 [42]. Ammonia, a chemical base, 
readily reacts with a number of acidic substances 
abundantly present in the atmosphere, frequently forming 
aerosols. It is also extremely soluble in water. 
Therefore, it is not surprising that fallout, washout, and 
rainout processes rapidly scavenge ammonia from the ait. 
Goncerrj for the amount of ammonia released into the 
atmosphere as a result of man's intervention into the 
nitrogen cycle would seem unjustified in light of such 
efficient mechanisms for removal of this gas. 
Unfortunately, some of the chemical products which result 
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in removal of gaseous ammonia from the atmosphere have 
deleterious effects on the environment. 

Acid Rain 

Atmospheric ammonia plays a key role in the chemical 
processes which lead to the formation of acid 
precipitation. Acid precipitation, hereafter referred to 
as acid rain, is defined as precipitation containing an 
excessive concentration of hydrogen ions. In terms of pH, 
acid rain is described as precipitation having pH below 
5.6, which is considered the minimal pH for pure water in 
the presence of atmospheric carbon dioxide [28] . A 
remarkable decline in the average pH of precipitation in 
North America and in Europe since the industrial revolution 
has been noted by numerous authors [10] US]. The most 
significant changes in pH have occurred within the last 
thirty years. The pH of rain and snow in many parts of the 
eastern United States and northern Europe now averages 
around 4.0; pH values as low as 2.1 have been measured in 
some storms [28]. Anthropogenically produced sulfur and 
nitrogen oxides, products of industry and combustion, are 
largely responsible for the drastic increase in the acidity 
of precipitation over widespread areas [1]. 

Chemical analyses of rainwater indicate that ammonium 
ions are present in significant quantities [9]. 
Researchers have examined possible chemical interactions 
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between ammonium and other ionic species present in 
raindrops and have reached somewhat surprising conclusions. 
A strong chemical base, ammonium acts to increase pH in 
raindrops. This effectively neutralizes some of the 
acidity caused by sulfate and nitrate ions also present in 
rainwater. However, other processes in which the ammonium 
ions participate apparently enhance the acidity of 

raindrops. According to Scott and Hobbs [39], Junge and 
Ryan [22], and others, ammonium acts as a chemical promoter 
within raindrops, speeding conversion of dissolved sulfur 
dioxide into acid sulfate. By controlling the pH of the 
system, ammonium ions also increase the solubility of 
gaseous sulfur dioxide in cloud droplets [12]. Ammonium's 
role in enhancing the solubility of gaseous sulfur dioxic. 
and promoting its oxidation to sulfate apparently 
supercedes its role as a neutralizing agent. Gaseous 
ammonia, which readily dissolves in raindrops to form 
ammonium ions, is now considered an important precursor to 
the phenomenon known as acid rain. 

Other Chemical Interactions and Their Effects 

Ammonia is known to undergo both aqueous-phase and 
gas-phase reactions with atmospheric pollutants to form 
aerosols such as ammonium nitrate and ammonium sulfate [I]. 
These particulates may affect the earth's radiation 
balance; more importantly, they constitute a health hazard 


in high concentrations. The Donora, Pennsylvania episode 
in 1948 is an ex^itmple of the health hazards brought about 
by such "smog." 

Ammonia molecules may also migrate to the st^ratosphere 
and participate in photochemical interactions with ozone, 
resulting in depletion of this essential stratospheric 
component. Such destruction of ozone molecules is thought 
to be of secondary importance compared to catalytic 
destruction of ozone by nitric oxide. At this time it is 
suspected that atmospheric ammonia may be a source for 
nitric oxide [32], but definite conclusions have not been 


reached . 
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II. LITERATURE REVIEW 

Early Measurements o£ Atmoiiplu^r ic Ammonia 

The role of ammonia in acid precipitation was not 
defined until Junge [22], Scott and Hobbs [39], and other 
pioneers in atmospheric chemistry published their results 
in the nineteen fifties and sixties. Therefore it is not 
surprising that attempts to measure gaseous ammonia were 
few and far between until fairly recently. The lack of 
interest in gaseous ammonia was also due in part to the 
widely held notion that ammonia was chiefly fixed to dust 
particles and did not occur in gaseous form at significant 
concentrations [21]. Early measurements of total ammonia 
(gaseous ammonia plus ammonium) in the nineteenth century 
reflect curiousity about the cycling of nitrogen compounds 
in the atmosphere, but after the turn of the century 
research all but ceased until the nineteen fifties [21]. 
In 1953 a network of sampling stations was established in 
Scandinavia. These stations produced monthly averages of 
gaseous ammonia concentrations. According to Junge [21] 
the results indicated a rather uniform distribution of 
ammonia — in the neighborhood of several micrograms per 
cubic meter (yg/m^)— with evidence of summer maxima and 
winter minima in ammonia levels at most stations. Junge 
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121J also conducted his own measurements during 1954 in 
Florida^ Hawaii, and Massachusetts. At these locations he 
obtained average ground-level concentrations of 5.1, 2.5, 
and 6.1 yg/m’, respectively. Junge tentatively concluded 
that the oceans are a source for gaseous ammonia, and 
suggested that thin organic films at the sea surface may 
release ammonia into the air. 

Almost twenty years later Georgii and Muller [11] 
published v,he results of their investigation of gaseous 
ammonia and ammonium-containing aerosols in the air over 
northern Europe. This study was the first comprehensive 
study of atmospheric ammonia. Prom 1969 until 1972 
monitoring was conducted over both land and water areas 
from aircraft outfitted with automatic sampling apparatus. 
Georgii and Muller reported that gaseous ammonia 
concentrations at ground level averaged 7 yg/m^ on "cold" 
days (temperatures below 10°C) and 18 yg/m^ on "warm" days 
(above 18°C). Seasonal differences in ammonia levels were 
still noticeable aloft.-^ At 3000 meters concentrations 
were, on the average, approximately 2 yg/m^ on cold days 
and 5 yg/m^ on warm days. Georgii and Muller found that 
ammonia concentrations generally reached constant 
"background" levels at 1500 meters on cold days and 3000 
meters on warm days. They concluded that atmospheric 
ammonia originated at the ground, and that temperature 
inversions as well as convection strongly affect the 
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vertical transport of ammonia. 

While most of the aircraft flights were centered over 
western Germany > measurements were also taken over the 
North Sea and over Bavaria. Ammonia levels were much lower 
over open water— as low as 0.3 pg/m Hence they concluded 
that ammonia is largely of continental origin. Georgii and 
Muller witnessed a classic example of rainout of ammonia 
over Bavaria. Upwind of their sampling area the 
intensification of a foehn wind system caused rainclouds to 
form over the Alps. Prom one day to the next ammonia 
concentrations in the lower troposphere fell from 20 yg/m ^ 
to 4 yg/m®. 

Georgii and Muller showed that meteorological 
parameters such as air temperature, lapse rates, and 
rainfall exert a strong influence on gaseous ammonia 
concentrations. No other investigation into atmospheric 
ammonia has been so thorough and informative. A number of 
other scientists, however, have made valuable contributions 
to our knowledge of atmospheric ammonia. For example, in 
1974 Lodge and coworkers [29] published results of their 
investigation into trace atmospheric constituents in the 
American tropics. Measurements were taken in 1967 and 
1968; these included seasonal and diurnal profiles. Lodge 
et al» found that gaseous ammonia concentrations averaged 15 
parts per billion (10 ) by volume (1.0 ppb=0.7 yg/m^ for 
ammonia at STP) , a level considerably higher than in 
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temperate climates. Their data indicated that ammonia 
concentrations were generally higher at night than during 
the day, perhaps reflecting inversion conditions. A 
"burst" of high concentrations at the beginning of the dry 
season was also noted by Lodge et al . 

Another important investigation was conducted by 
Tsunogai [41] in 1971. He found that the concentration of 
atmospheric ammonia was much lower over the 
ocean — averaging around 0.85 yg/m^ — than over land, where 
concentrations increased to approximately 3.4 y g/m^ . 
Tsunogai also measured the ratio of particulate aimnoniuro to 
gaseous ammonia over open ocean and compared it to data 
obtained near land. Prom his results he concluded that 
gaseous ammonia is mainly of continental origin and has an 
average residence time in the atmosphere of five to ten 
days [42] . 

Healy [15] [16] conducted several investigations of 
ammonia levels at both rural and urban locations in 
England. At rural Harwell he found that ammonia was 
present typically at 0.85-1.70 yg/m^. Healy concluded that 
domestic animals were responsible for most of the gaseous 
ammonia found over England. 

In the United States early measurements of gaseous 
ammonia were obtained by Breeding and his coworkers [6]. 
In 1971 and 1972 they conducted an investigation into 
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background levels of various trace atmospheric 
constituents. The sampling was done in rural areas around 
St. Louis, and attempts were made to pinpoint local 
sources of contamination. Breeding £t al. reported that 
background concentrations of ammonia varied between 2 and 6 
ppb in the area, with variations in this range due to 
"natural mechanisms." They detected no diurnal trends in 
ammonia levels. Interestingly, during a rainfall they 
measured higher concentrations than most of the samples 
obtained under sunny skies. 

Most early research into atmospheric ammonia in the 
United States was conducted in urban areas near strong 
sources of the gas [1). As would be expected, widely 
varying concentrations were reported. For example, one 
study obtained measurements as high as 450 ppb in the 
vicinity of a major dairy farm. Typically, it was reported 
that urban ammonia concentrations were significantly higher 
than those in rural areas. 

No federal or state agency in America has routinely 
monitored ambient ammonia. The unavailability of a 
continuous and reliable method for measuring ammonia at low 
concentrations contributed to the paucity of measurements 
[11 . 
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Recent Measurements Using New Techniques 

Continued interest in gaseous ammonia has spurred 
development of more sensitive and accurate measurement 
techniques. Most of these new methods are still in 
experimental development, but very recently results of some 
of these measurements have begun to appear in scientific 
publications. Many of these reports involve measurements 
taken in Southeastern Virg inia--both by remote and situ 
techniques--and thus form an excellent basis for comparison 
with the measurements which will be reported in this paper. 

Using a new ^ situ technique in 1978, McClenney and 
Bennett [30] measured gaseous ammonia at two sites in North 
Carolina. Near Research Triangle Park, where five 
measurements were made, the average ammonia concentration 
was 4.1 ppb. McClenney and Bennett also sampled at Cedar 
Island, North Carolina. No attempt was made to delineate 
diurnal trends in ammonia levels (all sampling took place 
in the afternoon) , but an effort was made to sample near 
possible local sources of ammonia. Samples were taken from 
a ten foot height over land (average ammonia concentration: 
1.0 ppb), from a height of six feet over marsh grass (one 
sample: 0.9 ppb), from a height of two inches over marsh 
grass at the water's edge (average: 2.9 ppb), over land at 
a height of two feet (0.3 ppb), and from within a shelter 
during rain (average concentration: 0.3 for six samples). 
McClenney and Bennett avoided drawing any broad conclusions 


Page 14 


from such spares data. Moreover, their main objective was 
to test the efficacy of their sampling methods and 
analytical technique. , 

A team of reseaTchers at NASA Langley in Hampton, 
Virginia and at Old Dominion University began conducting 
both ^ situ and remote measurements of ambient ammonia in 
1979. One of their reports [18] compared ^ situ data 
obtained from an aircraft measurement platform with 
vertical profiles of ammonia concentrations from a Remote 
Infrared Heterodyne Radiometer (IHR). Results from the in 
situ data were in good agreement with the remotely obtained 
profiles. Both techniques indicated that on the day the 
measurements took place (June 12 , 1979) ammonia levels were 
approximately 1 ppb near the ground and slowly decreased 
with increasing altitude. 

In another report (19J these researchers compared 
vertical ammonia profiles taken in March of 1979 with 
profiles obtained in August of that year. They found that 
ammonia concentrations were significantly lower throughout 
the lower troposphere in August than in March. The latter 
measurement indicated that ground-level concentrations of 
gaseous ammonia were approximately 1.3 ppb, while March 
levels were around 10 ppb. These results were contrary to 
Georg ii and Muller's in situ measurements in Europe, which 
indicated that ammonia release from the soil was dependent 
on soil temperature. The measurements, however. 




qualitatively agreed with measurements obtained on Long 
Island in 1976 by Peyton et al. [36] during a preliminary 
test of a similar IHR instrument. Peyton and coworkers had 
found that ammonia levels near the ground had decreased 
from 12 ppb in March to 0.2 ppb in July. These results 
indicated to the researchers at Langley that other factors 
beside temperature influenced ammonia emission from the 
soil. It was concluded that rapid volatilization of 
ammonium nitrate fertilizer which had been applied to 
nearby fields several weeks earlier had contributed to the 
high ammonia concentrations found at Langley in March of 
1979 [271. 

Extensive ^ situ measurements in 1979 and 1980 at 
Langley (some of which overlapped this study's research) 
indicated to Hoell £t £l. [20] that soil moisture was a 
factor which also influenced local ammonia levels in the 
>?^tmosphere . They found that background ammonia levels in 
the late summer of 1980 were much lower than during the 
corresponding period of 1979. Abundant rain fell in both 
August and September of 1979, when ammonia levels averaged 
nearly 2 ppb. In 1980, hqis'ever, severe drought apparently 
contributed to much lower ammonia concentrations, which 
averaged only around 0.5 ppb in August and September. It 
is also noteworthy that the spring maximum in 1980 was much 
less pronounced with the in situ data than would be 


expected from the previously reported March 1979 levels 
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recorded by remote means. The authors noted that the 
entire year of 1980 was relatively dry^ and hence a variety 
of climatological conditions may be operating 
simultaneously to influence concentrations of ammonia in 
the atmosphere. 
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III. STATEMENT OF THE PROBLEM 


Current Requirements for ^ Situ Ammonia Measurements 

Research into the distribution and chemical 
interactions of atmospheric ammonia has been hampered by a 
lack of data concerning this important trace gas. Most of 
what we know about ammonia in the atmosphere comes from the 
handful of articles and reports outlined in the foregoing 
discussion. In the past few years federal regulatory 
agencies, atmospheric chemists, and other concerned groups 
and individuals have expressed the desire for more 
information on the temporal and spatial distribution of 
gaseous ammonia, its sources and sinks, its residence time 
in the atmosphere, and concentrations of related trace 
species. For example, in 1979 the Subcommittee on Ammonia 
(of the Committee on Medical and Biologic Effects on 
Environmental Pollutants, National Research Council) 
recommended that accurate estimates of the emission, 
movement, and degradation of ammonia in the atmosphere be 
obtained. It also listed as a high priority further study 
of ammonia's complicated atmospheric chemistry. Toward 
this end: 
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Simultaneous measurement of ammonia and of particulate 
hydrogen (acidity), ammonium, sulfate, and nitrate content 
are needed to elucidate further the role of ammonia in the 
formation of |>articulate ammonium, nitrate, and sulfate and 
to formulate improved strategies for the control of these 
major inorganic pollutants (1). 


Furthei.nore, the Subcommittee recognized that current 
monitoring methods for ammonia are woefully inadeqj ite and 
stated that "methods should be developed or refined for the 
routine measurement of ambient ammonia at parts-per-bill ion 
concentrations. These methods should be suitable for 
continuous measurement of ambient ammonia as part of a 
limited monitoring network." [IJ 

The existence of a monitoring network would be a boon 
to scientists interested in modeling the complicated 
atmospheric chemistry of ammonia. Under non-steady-state 
conditions, validation of a model's output is impossible 
without extensive spatially and temporally resolved data. 
Under the simplifying (if somewhat unrealistic) assumption 
of steady-state, horizontally homogeneous distributions of 
gaseous ammonia, modelers can "make do" with vertical 
profiles of ammonia concentrations. In fact, vertical 
profiles obtained by IHR have been used by a team at NASA 
Langley in conjunction with a steady-state model they 
developed. Using the data thay were able to derive rate 
constants for heterogeneous reactions (including gas to 
particle conversions, dry deposition, and rainout) 
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involving ammonia [27]. 

At the present time, the IHR can yield data useful to 
numerical modelers at Langley with less cost and bother 
than available ^ situ techniques. However, jUn situ 
monitors are still needed to provide ground truthing of the 
data f om the THR and to extend measurement capabilities. 
The IHR system at Langley has limited sensitivity (0.5 ppb) 
and cannot be used when direct sunlight is not present 
[ 20 ] . 

An automated Jji situ ammonia monitor is the only 
practical method for obtaining continuous data. If such a 
monitor were to exist, it could be easily integrated into 
existing monitoring networks which have been established by 
various state and federal agencies. Without such a network 
scientists cannot accurately pinpoint sources of the gas or 
reach definite conclusions about the influence of various 
meteorological parameters on the temporal and spatial 
distribution of ammonia and related atmospheric species. 
While the utility of a single, isolated monitoring station 
is limited, inferences can be drawn where conclusions 
cannot. The present data base is so sparce that any 
information about ammonia, however fragmentary, would be 
welcomed by the scientific community. 
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Available Measurement Techniques 

Most methods of measuring atmospheric ammonia involve 
a period of sampling (by bubbling or some other form of 
preconcentration), followed by an analysis of ammonia 
content in the sample. Collection of ammonia is 
complicated by a number of difficulties. Contamination of 
samples by ammonia emanated from nearby personnel is a 
problem which many investigators have encountered. Because 
ammonia adsorbs onto almost any surface, care must also be 
taken to minimize contact of the air stream containing 
ammonia with the walls of tubing. An additional problem 
arises out of the fact that ammonia and ammonium aerosols 
are generally both present in air. Differentiation between 
the two is impossible with most analytical techniques. The 
inclusion of particulate prefilters in the sampling stream 
may result in either positive or negative errors in 
measurement, depending on the aerosol content of the air 
and the type of filter used [1|. Sampling efficiency is 
yet another problem. Most sampling techniques rely on the 
high solubility of ammonia in acidified water. But studies 
have shown that bubblers are usually less than 100% 
efficient in dissolving ammonia. Efficiency also decreases 
with low ammonia concentrations [1]. 

Assuming valid samples of ammonia have been obtained, 
any of several wet chemical techniques can be used to 
analyze the samples for ammonia. The Nessler method. 
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involving the reaction of ammonium ion with Nessler's 
reagent, is considered the classic colorimetric technique. 
It has been demonstrated, however, that interfering 
substances frequently invalidate the measurement of ammonia 
by this method [34]. Alternatives to the Nessler method 
have been developed and widely used in recent years. The 
two most popular colorimetric techniques are the indophenol 
method and the pyridine-pyrazolone method. Breeding et £l. 
[6] and Georgii and Muller [11] relied or* the former 
technique to determine ammonia concentrations for their 
extensive investigations. Okita and Kanamori [34] have 
extolled the virtues of the latter method. 

Several reservations about these and other aqueous 
techniques immediately come to mind. As mentioned earlier, 
none of the aqueous methods is capable of separating the 
gaseous ammonia from species containing ammonium ions. 
Indeed, to avoid the problems inherent in the use of 
prefilters, some investigators have contented themselves 
with simply measuring the total ammonium concentrations 
present in the air. The three techniques discussed above 
are generally considered among the most sensitive of all 
colorimetric methods. Yet their sensitivity limits lie in 
the range of 0.01-0.05 mg/liter of solution [1]. Using 
standard bubbling apparatus, less than ten liters per 
minute (1/min) of air per liter of solution can be sampled 
with high efficiency. This implies sampling times of two 
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hours or longer for valid determination of background 
ammonia concentrations (less than ten nanograms/liter). 
Breeding et ^1, (6) sampled fo,r one hour, but they gave 
standard deviations for accuracy of the measurements larger 
than the approximate average concfintrat ion of ammonia in 
many cases. Other researchers report that a sampling 
period of anywhere between one and two hours was necessary 
to obtain measureable amounts of ammonia [11] [34] . 

Other wet chemical analytical methods reported in the 
literature include the specific-ion electrode method and 
ion chromatography [1]. Both claim sensitivity limits only 
slightly greater than the colorimetric techniques 
(approximately 0.1 mg/liter of solution). Moreover, 
bubbling apparatus must still be used for sampling before 
employing either of these means of analysis. 

To avoid the pitfalls associated with bubbling 
apparatus, Shendrikar and Lodge (40) adapted the so-called 
ring oven technique for analysis of ammonia. Basically the 
sampling technique involves the impinging of ammonia (and 
ammonium ions as well) on a filter paper impregnated with 
oxalic acid. Pinal determination of ammonia content 
depends on visual comparisons of sample spots on the paper 
with the intensities of spots obtained from known amounts 
of ammonia. Not only is the procedure complex, but long 
sampling times (approximately two hours) are necessary to 
obtain measureable quantities of ammonia. Interference by 
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formaldehyde is also a problem, according to the authors. 

It is readily apparent that all of the wet and 
quasi-wet analytical techniques have their shortcomings. 
The most fundamental of these which all of the above share 
is that ammonia in its gas phase is never measured 
directly, but is instead converted into a form (ammonium 
ion) which may be detected by chemical means. Many direct 
methods of measuring gaseous ammonia now exist, but almost 
without exception these lack the sensitivity to accurately 
assess the minute quantities of ammonia present in the air 
of most environments. 

One possible exception was alluded to earlier. The 
technique involves the use of an IHR, which measures the 
solar absorption due to atmospheric ammonia. One such 
system is being used at the NASA Langley Research Center. 
Its current sensitivity limit (0.5 ppb) makes it only 
marginally suitable for detecting background levels of 
ammonia. However, improvements in sensitivity are possible 
[13]. Since it relies on the sun for a radiation source, 
the IHR's use is restricted to relatively clear days. 
Offsetting this disadvantage is the system's capability to 
obtain virtually instantaneous profiles of ammonia 
concentrations from ground level up to 30 km [19] . 

Among the other direct techniques which have been 
experimentally employed in attempts to measure ambient 
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ammonia are Four ier-transforra long-path infrared 
spectroscopy, second derivative spectroscopy, 
chemiluminescence, and laser photo-acoustic spectroscopy 
(LPS) (11. These techniques currently cannot detect low 
background concentrations of ammonia. However, new methods 
of preconcentrating ambient ammonia now exist which make 
routine monitoring of atmospheric ammonia by either 
chemiluminescence or LPS practical. It was stated earlier 
that gaseous ammonia has the propensity to adhere to almost 
any available surface. Ironically, this characteristic has 
opened up new avenues of research into substances which 
adsorb and subsequently can be induced to desorb ammonia 
molecules . 

Two substances, tungsten oxide and small teflon beads 
(known by the trade name "Chromosorb T" ) , have thus far 
been found which are suitable for preconcentrating ammonia. 
Both substances exhibit a strong tendency to "trap" ammonia 
molecules on their surfaces at room temperature; upon 
heating they release the ammonia at an accelerated rate. 
Tungsten oxide has emerged as the preferred substance for 
most purposes. The teflon beads are easily fused by 
excessive temperatures, so the ammonia may not be driven 
off as rapidly as is the case with tungsten oxide. This 
results in a lower factor of preconcentration for the 
teflon beads. Furthermore, the teflon beads are packed in 
quartz tubing while the tungsten oxide merely coats the 
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Inside of the tubes. Ammonium-containing particles pass 
through the tungsten oxide-coated tubes so long as the flow 
is nearly laminar within the tubes [5]. The particles must 
be prefiltered from the teflon bead-packed tubes, and this 
introduces the uncertainties alluded to earlier. 

The ammonia, once driven off from the tubes, may be 
analyzed by a specially modified oxides of nitrogen 
analyzer, which utilyzes the chemiluminescent reaction 
between nitric oxide and ozone. The ammonia molecules are 
first converted to nitric oxide by a catalyst (platinum at 
1000 °C is most frequently employed) before being passed 
through the analyzer. Unfortunately, many other nitrogen 
compounds which may be present in the air stream are also 
converted. Studies have shown that nitric acid and some 
amines are collected by both the teflon-packed and tungsten 
oxide-coated tubes, and these constitute a positive 
interference when the tubes are used with a 
chemiluminescent oxides of nitrogen analyzer (5). 

The chemiluminescent technique uses commercially 
available gas analyzers which have been specially modified. 
These instruments are easy to use and relatively portable. 
In addition, excellent sensitivity has been reported [5). 
These advantages are offset by the lack of selectivity of 
the chemiluminescent oxides of nitrogen analyzers, which 
allows possible interference by a number of common 
atmospheric contaminants. 
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A more specific physical method for determination of 
ammonia is still in experimental development. A brief 
history of the discovery and development of the principles 
underlying this technique follows. 

Laser Photo- acoustic Spectroscopic Measurement of 
Atmospheric Ammonia 

The science of acoustics owes a great deal to the work 
of Bell, Tyndall, and Roentgen in the late nineteenth 
century. In 1881 they first reported on a phenomenon which 
became known as the photo-acoustic effect (23]. As its 
name implies, the photo-acoustic effect occurs when some of 
the energy absorbed by gas molecules from a beam of 
radiation results in the net heating of the gas. If the 
beam is made to vary in intensity, the temperature of the 
gas in turn fluctuates. As a result the pressure exerted 
by the gas will undergo fluctuations if the volume of the 
gas is held constant. Pressure fluctuations in the gas 
medium propogate as sound waves. A pressure transducer, 
which is no more than a sensitive microphone, may be used 
to convert the sound into electrical signals. 

Many uses have been found for the photo-acoustic 
effect since its discovery. Commercial detectors which use 
the photo-acoustic effect to measure carbon dioxide, carbon 
monoxide, or water vapor have been around for a number of 
years [23] . They rely on conventional black-body emitters 
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of radiation and employ optical filters to select the 
region of the electro-magnetic spectrum where the gas being 
measured strongly absorbs radiation. Ammonia is generally 
too scarce to be successfully monitored by such methods, 
although it does strongly absorb certain wavelengths of 
infrared radiation. Other gas species which absorb in the 
nearby wavelengths effectively "drown out" any slight 
response which may be due to ammonia molecules alone. But 
gas lasers produce virtually monochromatic radiation. 
Isolated absorption features (single absorption lines 
rather than wide bands in which absorption takes place) can 


be used to 

detect minute 

amounts of 

absorbing gas. 

The 

technique 

employing 

gas lasers 

to produce 

the 

photo-acoustic effect 

has come to 

be known as 

laser 


photo-acoustic spectroscopy. Its rapid development in the 
last ten years is attributable to the increasingly powerful 
and reliable carbon dioxide lasers which have only recently 
become commercially available. 

One of the pioneers in the application of laser 
photo-acoustic spectroscopy (LPS) to pollutant monitoring 
has been L.B. Kreuzer. In the early nineteen seventies he 
developed a prototype multiple gas analyzer which relied on 
a discretely tunable carbon dioxide gas laser and a 
photo-acoustic effect detector cell [23]. The gas laser 
could be tuned to a number of wavelengths which roughly 
correspond to absorption lines of pollutant molecules. 


Page 28 


Multiple pollutant gas analysis was accomplished by solving 
a set o£ simultaneous linear equations containing 
information on the signal response, absorptivity of 
individual gases, path length, and power of the infrared 
beam at each wavelength. This chore was done by a 
microcomputer. Uncertainties in the calculation of 
absorptivities and limitations to the sensitivity of the 
detector cell limited the instrument's aensitivity to the 1 
ppb range for most gases. Still, for ammonia at least, 
this was a considerable improvement over previously 
existing systems. 

Kreuzer's shotgun approach to pollutant detection (the 
laser was tuned to 200 different wavelengths in a five 
minute time span) was aimed at replacing a whole bank of 
gas analyzers with one device which could almost 
simultaneously measure up to ten gases. A team at NASA 
Langley working with an EPA group from Research Triangle 
Park, North Carolina needed a more sensitive instrument for 
measuring background concentrations of gaseous ammonia. 
They developed their own LPS system in the late nineteen 
seventies for this purpose. The use of preconcentrator 
tubes and other technological refinements incorporated into 
this LPS system make it possible to monitor ammonia 
concentrations of less than 0.1 ppb on an hourly basis. 
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IV. PRINCIPAL FOCUS OF THIS STUDY 

Parameters of Interest Relating to the Measurements 

A measurement program for any trace atmospheric 
constituent gains significance if attempts are made to 
relate the measurements to parameters which may have a 
bearing on its distribution in space and time. An attempt 
will be made in this paper to analyze ammonia 
concentrations in terms of several meteorological, 
climatological and soil parameters. Limitations in the 
data may prevent positive conclusions concerning possible 
cause-and-ef feet relationships, but inferences can 
certainly be drawn where conclusions cannot. 

Therefore the ammonia measurements will be categorized 
by ambient temperature intervals, in order that a possible 
relationship between ammonia concentrations and air 
temperature may be investigated. In addition, a regression 
analysis of ammonia concentrations versus air temperature 
will be performed. 

Further, the data will be divided into eight groups 
corresponding to the prevailing wind directions at the 
times of sampling. In this way the suspected continental 
origin of most gaseous ammonia may be demonstrated. 
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Perhaps local sources of the gas can be inferred from the 
data as well. 

Seasonal variations in atnmonia levels will be 
described, the results will be compared to other ammonia 
measurements obtained at Langley, such as those reported by 
Hoell et al. (20]. If it is determined that local ammonia 
levels repeatedly peak in the early spring, as Hoell £t al. 
claim, a possible explanation for the peculiar early 
springtime maximum will be further investigated. 

Other data obtained by previous measurement programs 
at NASA Langley will be compared to data from this program 
to ascertain if systematic variations exist in ammonia 
levels from year to year. If such discrepancies do occur, 
they may be due to differences in soil moisture. A soil 
water budget model has been developed at Old Dominion 
University which will be used to estimate soil moisture 
levels in both 1979 and 1980, In this manner a causitive 
link between moisture levels in nearby soils and local 
ammonia production may be inferred. 

Moisture in the air may also affect ammonia 
concentrations. The measurements will be analyzed in terms 
of both water vapor mixing ratio and relative humidity in 
an attempt to ascertain what, if any, relationship exists 
between ammonia concentrations and moisture content of the 


air. 
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An evaluation of the Merits and Limitations of This 
Technique 

Several promising new methods for measuring 
atiAospheric ammonia have emerged in recent years, but none 
of these has yet been proven superior for all purposes or 
situations. Remote techniques such as infrared heterodyne 
radiometry, while they hold great promise, may never be 
able to totally supplant ^ situ methods, especially where 
continuous monitoring is desired. Modified 
chemiluminescent oxides of nitrogen monitors may now 
possess the necessary sensitivity and selectivity to 
differentiate gaseous ammonia from other ambient nitrogen 
compounds, but continuous, routine measurement using 
chemiluminescence has not yet become a reality. The same 
is true of laser photo-acoustic spectroscopy-. The method 
is still in experimental development and further 
improvements are needed. A second focus of this paper will 
therefore center on an evaluation of the system’s 
performance during the measurement program. The aim will 
be to determine where further improvements ace necessary 
and to suggest possible modifications to the existing 
design or measurement procedure. Toward this end a 
narrative discussion of the problems and discoveries 
encountered during the course of the research will be 
presented. 
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V. THE MEASUREMENT SYSTEM 
Theoretical Background 

The ammonia molecule has a pyramidal structure and 
behaves as a symmetr ic-top rotor [1). Whi],e the gas is 
virtually transparent in the visible and near-ultraviolet 
regions of the electromagnetic spectrum, it exhibits strong 
absorption in the near-infrared. In the ten micron 
wavelength region particularly, the gas has large 
absorption coefficients 1,371. This spectral region 
corresponds to a series of strong absorption bands caused 
by the molecule's vibrational-rotational characteristics. 
Figure 1 shows line strengths of individual lines for 
ammonia and several other molecules. Several strong 
absorption lines for ammonia are present in the region 
around 10.78 microns, or 927 Kaysers (1 Kayser » 1 cm-1). 
In order to isolate absorption of radiation due to the 
presence of ammonia from absorption by other molecules 
which may be present, a strong absorption feature which 
does not significantly overlap with the absorption lines of 
other molecules * must be isolated. One such absorption 
feature is centered at 92 7. 32323 (+0. 0001 2 ) Kaysers. This 
corresponds to the aQ(6,6) 2 fundamental 
vibrational-rotational transition of the ammonia molecule 
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(17J. The absorption feature has a line strength of 
4.09*10-19 Kaysers/( molecule cm**2) [19]. Nearby 
absorption lines of other gases are generally several 
orders of magnitude weaker. Water vapor, a much more 
abundant gas in the atmosphere, may cause interference when 
ambient air is being analyzed for ammonia. The sampling 
method used in this study, however, virtually eliminates 
the possibility of water vapor causing errors. Evidence 
for this assertion will be presented in another section. 

Since the spectrum is crowded with absorption lines, a 
means must be found to provide radiation across a very 
narrow spectral interval. Lasers provide virtually 
monochromatic radiation of great intensity. For this LPS 
system an isotopic carbon dioxide laser (carbon 13, uxygen 
16) was chosen which lases continuously with power on the 
order of several watts at S27. 300406 Kaysers [19]. (This 
corresponds to the Rl 8 line in the (00° 1-10° 0) band.) Since 
the half-width of the absorption feature of interest is 
0.08 Kaysers [19], and the incident radiation is within 
0.023 Kaysers of the line center, strong absorption is 
assurtd. 

When dealing with minute amounts of ammonia, direct 
measurement of absorption (or beam attenuation) is 
impractical, if not impossible. Fortunately an indirect 
means of determining absorption exists: the photo-acoustic 
effect. Once incident radiation has been absorbed by an 
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ammonia molecule, there is a large probability that 
collisional deactivation will occur before the molecule has 
a chance to reradiate the absorbed energy [7). Collisional 
deactivation occurs at normal atmospheric pressures and 
temperatures because of the large number of molecular 
collisions taking place over a short time in the crowded 
gas medium. Excited molecules donate their increased 
vibrational-rotational energy to the colliding molecules, 
resulting in larger kinetic energies for the latter. When 
the incoming radiation is "chopped" (made to fluctuate), 
thermal energy and pressure within the detector cell in 
turn fluctuate. Sound waves are thus produced which 
propagate through the medium. 


For small concentrations of ammonia , pressure waves 
resulting from absorption of radiation could easily be 
obscured by the overall "noise" within the detector cell. 
Additional steps must be taken to enhance and isolate the 
sound produced by the action of the ammonia molecules 
alone. The photo-acoustic cell may be so designed to 
operate as a Helmholtz resonator at certain sound 
frequencies. The addition of a lock-in amplifier to the 
system further enhances the sensitivity by reducing the 
bandwidth . 


A Helmholtz resonator is an acoustic enclosure whose 
geometry allows the gas within a cavity to "communicate" 
with a second chamber through a narrow neck. Viscous 
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dissipation of the sound waves in the secondary chamber is 
thus minimized because they "feed back" into the main 
cavity of the resonator, which returns most of the energy 
to the neck in the form of reflected waves. There is a 
fundamental resonant frequency for a given resonator, which 
is a function of the speed of sound, the cross sectional 
area and length of the neck, and the volume and geometry of 
the cavities (38]. In the LPS oystem the light chopper at 
the proper frequency creates resonance within the detector 
cell . 


The lock-in amplifier is a specialized ac amplifier 
which synchronously demodulates only certain frequencies of 
an incoming "noisy" signal at a reference frequency. It 
filters out electrical impulses not in phase with or at the 
same frequency as the reference signal. An electrical 
signal from the mechanical light chopper acts as a 
reference signal. Since the microphone signals have a 
frequency component (due to molecular absorption) that is 
the same as that of the light chopper, the signal-to-noise 
ratio of the detection is greatly increased. 

Theoretical limits to the amount of ammonia which may 
be detected by LPS are eventually encountered and must bt 
mentioned briefly here. After Kreuzer [23], the signal (S) 
which is measured by the microphone within the detector may 
be mathematically expressed as a function of the laser 
power (P) and the absorbance (A) of the ammonia molecules; 
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S»P* (l-exp(-A) ] . This expression shows that the signal is 
directly proportional to the laser power. The exponential 
term results from Peers 's Law, The absorbance in turn 
effects the power transmitted through the sample (T): 

' T“P*exp(-A) . Combining these equations and solving for 
sample absorbance; A»ln(S/T+l); and, according to Beer's 
Law, the absorbance may also be described as a function of 
the path length 1, absorptivity of the gas a, and gas 
concentration C: A=laC, Therefore, the concentration of 

ammonia is given by C»( a/1) *ln(S/T+l ) . 

This last equation defines the basic parameters which 
must be dealt with in designing a LPS system. Moreover, 
since the detector used in this study operates as a 
Helmholtz resonator, other factors also influence the 
sensitivity of the apparatus. These include the geometric 
configuration of the resonator, connecting neck, and 
absorption chamber, the frequency of the chopped radiation, 
the molecular viscosity coefficient of the gas, and the 
speed of sound within the gas medium [31]. The combination 
of photo'-acoustic cell and Helmholtz resonator increases 
the system's responsivity to absorption by ammonia 
molecules, but the "window signal" due to absorption of 
radiation at interior surfaces of the cell is also 
amplified. Therefore the magnitude of the "window signal" 
essentially establishes the minimum detectable limit for 
this system [31]. Currently around four nanograms of 
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ammonia can 

be 

detected by 

the LPS 

system. 

This is 

adequate 

for 

obtaining 

hourly 

averaged 

ammonia 


concentrations under almost any conditions. 

Apparatus 

In its present stage of development the LPS system 
used for this study does not directly sample ambient air 
for determination of ammonia concentration. Measurements 
are made in a two step process; a period of sampling or 
preconcentration of ammonia is followed by subsequent 
analysis of the sample for ammonia content. It is 
therefore convenient to describe the sampling equipment, 
analytical apparatus, and calibration system separately. 

Limitations in the sensitivity of the analytical 
apparatus necessitate some method of preconcentration of 
the ambient ammonia before analysis. Hollow quartz 
collection tubes coated on the inside with tungsten oxide 
have been used exclusively during the course of this study. 
The tubes are approximately 40 cm long, 35 cm of which has 
been coated. The inside diameters measure approximately 
0.4 cm. These tubes were developed by Dr. Robert Braman 
of the University of South Florida during a similar 
measurement program. A likely mechanism for the action of 
the tubes in selectively adsorbing ammonia has been 
described by Braman (5J. The tungsten oxide interacts with 
water vapor to create tungstic acid (H^WOj^). 


This weak 
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acid in turn reacts with gaseous ammonia, forming NH4HW04> 
At temperatures above 350°C the reaction is reversed, 
liberating the ammonia. The tubes are wrapped with 
nichrome wire so that electrical power may be applied to 
heat the tubes at the time of analysis. More observations 
on the tubes' performance will be presented in a later 
section. 

The remainder of the sampling apparatus consists of a 
suction pump and associated tubing with three calibrated 
flowmeters and control valves incorporated into the air 
stream. Multiple valves and flowmeters allow replicates to 
be taken. The entire sampling apparatus is portable and 
requires only standard 110 Vac electrical outlets for 
operation. Sampling times are monitored by means of a 
stopwatch. 

If the sampling system is simple and compact, the 
analytical apparatus compensates for this by virtue of its 
complexity and bulkiness. Figure 2 is a simplified diagram 
of the basic components of the LPS analytical system. 
Mechanically chopped (approximately 1000 Hz) radiation from 
the laser (GTE Sylvania Model 950) is directed by a series 
of mirrors and a focusing lens through the narrow opening 
of the photo-acoustic detector. After passing through the 
detector, the beam is divided by a beam splitter. A 
portion of the beam is directed into a power meter, as 
shown in the diagram. The remaining radiation finds its 
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ZERO AIR EXIT AIR 



Figure 2. Simplified diagram of the basic components of 
the LPS analytical system. 
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way into a spectrum analyzer (Optical Engineering C02 
Spectrum Analyzer), so that the wavelength of laser 
radiation may be monitored and adjusted, if necessary. 
Ammonia passes into the detection cell directly from the 
heated tubes. Helium (technical grade: 99.995% pure) is 
used as a carrier gas (flowrate: 0.1 1/min) because of its 
inert and nonabsorbing qualities. The photo-acoustic 
detector is coated with teflon and is heated (approximately 
50°C) to minimize retention of ammonia by its walls. 
Dimensions of the cell through which the gas stream passes 
are approximately 25 cm long by 0.6 cm in diameter. The 
microphone (Bruel and Kjoer Model 4138 with Type 2801 power 
supply) is mounted on top of the Helmholtz resonator. Its 
signal is synchronously demodulated by the lock-in 
amplifier (EG and G Princeton Applied Research Model 510) 
with the chopper frequency as a reference. Any signal in 
phase with the chopper is amplified and fed to the 
recorder. Simultaneously the laser power as measured by 
the power meter (Molelectron Corp. Pyroelectric 
Radiometer) is also recorded. 

Early in this study it became apparent that 
fluctuations in laser power produced signal fluctuations 
from the photo-acoustic detector, creating noise and 
producing large errors in the measurements. This problem 
was eventually solved by including a ratiometer (Princeton 
Applied Research Model 193 Multiplier/divider) in the 
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processing system. This instrument electronically divides 
the signal from the detector by the power reading from the 
power meter. Inclusion of the ratiometer greatly enhanced 
the performance of the system and allowed accurate 
measurement and calibration even when the laser power* was 
markedly oscillating. The "transparent" apertures at both 
ends of the photo-acoustic detector created a slight 
"window signal" because of absorption of radiation by the 
aperture material. This constituted no great problem, 
since the recorder could be offset to compensate for the 
window signal. 

The two channel recorder (Linear Instruments Corp. 
Integrator/recorder Model 282) is equipped with an 
electronic integrator with adjustable baseline response and 
three "count rate" settings. The count rate recorded by 
the electronic integrator is governed by the signal 
received from the lock-in amplifier. As ammonia is driven 
off from the collector tubes, a curve is described on the 
recorder paper. The area under this curve before it 
returns to baseline is proportional to the amount of 
ammonia which passed through the detector. The integrator 
counts the number of area units contained within the 
envelope of the curve, eliminating the need for determining 
the area by graphical means. 

The remainder of the system's components consists of 
calibration equipment. See figure 3 for a diagram of the 



Figure 3. Diagram of the calibration system. 
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calibration system. Sorbants are used to clean and dry 
pressurized "house" air# which then passes into a triple 
dil^ition system. Ammonia is added to the air by means of a 
commercial double dilution device (Metronics Model 340 
Oynacal ibrator ) containing an ammonia permeation tube 
(Metronics wafer type) housed in a temperature-controlled 
chamber. Triple dilution is achieved by the addition of a 
flowmeter and control valve to regulate another air stream. 
Clean house air is humidified while passing through this 
third circuit by inclusion of a bubbler containing 
distilled water. Humid air is used (relative humidity: 
60%) because this at least partially duplicates the ambient 
environment, and because the sampling tubes apparently 
require at least a minute amount of water vapor to maintain 
their sensitivity. (Water vapor must be available to allow 
conversion of tungsten oxide to tungstic acid.) 

Other calibration equipment includes a timer connected 
to a shutoff valve to control calibration sampling time, a 
pump used to pull air through the tubes when performing 
calibration sampling, a valve and flowmeter to regulate the 
flow through the tubes, and various connectors and brackets 
for mounting the tubes. Tubes are mounted in a teflon 
"tee" connector during calibration sampling. One outlet of 
the "tee" remains open to the room. This prevents a 
pressure buildup due to differences between the flowrates 
of the calibration air stream and the sampling air stream. 
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Since the calibration air Si^ream is over ten tiroes that 
demanded by the tubes when samping takes places no backwash 
of room air can take place through the open "tee." 

Cal ibration 

Calibration of the instrument was accomplished by 
linear regression analysis of the integrated recorder 
response to several different "loads" of ammonia in a 
collection tube. Each collection tube was separately 
calibrated, although frequently different tubes' regression 
equations were very similar. Absolute calibration was 
traced to the permeation of ammonia through teflon (33J . 
Since the permeation rate is dependent on air temperature, 
the permeation device was kept at a constant temperature of 
30°C throughout this study, except for brief periods when 
the device was being weighed. Periodic weighings allowed 
the permeation rate to be determined gravimetrically . 
Figure 4 shows the data, best linear fit, and regression 
equation for a permeation device used during most of this 
study. The nearly linear weight loss (linear correlation 
coefficient; -0.9999051) corresponds to an average elution 
rate of 58.664452 nanograms per minute (ng/min). For the 
great majority of calibration runs flow rates through the 
permeation chamber, the secondary dilution, and the 
tertiary (humidified) dilution were fixed at 0.333, 5.420, 
and 8.70 1/min, respectively. This resulted in an ammonia 
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concentration of 4.05 ng/1 in the calibration air stream. 
A propagation of errors analysis has been performed which 
estimates the magnitude of errors in the gas flow 
measurements. The analysis indicates that the 
concentration may have varied by +3.5% from the 4.05 ng/1 
mean value. For simplicity in calculations the 
concentration was assumed to be 4.0 ng/1. Ammonia loading 
was varied by pulling 1 1/min of the calibration mixture 
through the tubes for various periods of time. For 
example, a 24 ng load of ammonia was achieved by sampling 
for six minutes. Since the tubes tended to exhibit some 
changes) in response over varying lengths of time, frequent 
recalibrations were performed on the tubes to maintain 
accuracy. Frequently chocks were also made after 
measurement runs to ensure that sudden response changes had 
not occured. VJhenever it was felt that a new calibration 
was required for a given tube, at least five and often more 
runs at various loads were made in order to derive a new 
regression equation. Linearity of the regression equations 
was generally excellent; linear correlation coefficients 
averaged over 0.98 for some 60 calibrations. Tubes that 
exhibited poor reproducibility were not used for 
measurements until reproducibility was restored. Linear 
correlation coefficients less than 0.9 (95% confidence 
level for 5 data pairs) were considered unacceptable. 
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In order to arrive at an estimate of the errors in 
deriving regression equations^ a "worst case" experiment 
has been performed for an individual tube. Several 
individual calibrations were lumped together into an 
overall regression analysis covering several months of 
sampling. Figure 5 illustrates the results of this 
experiment. An envelope representing the slope and 
intercept errors (two standard deviations from the mean) 
has been drawn around the regression line. Inspection 
reveals the uncertainty in the slope and intercept values 
to be tolerably small (less than 2.5% uncertainty in the 
slope). However, for small loads tne relative errors in 
calibration might be proportionally larger than for large 
loads. Another experiment confirmed this suspicion. A 


calibrated tube 

that 

had 

been used for measurement 

was 

later subjected 

to 

a follow-up , analysis 

of 

its 

reproducibil ity 

for 

small 

loadings of ammonia. 

It 

was 


found that the mean of the area units corresponded to a 7.7 
ng average loading, according to the previous calibration. 
This systematic error implies a slight loss of sensitivity, 
although the standard deviation about this mean value — 0.6 
ng — was large enough to indicate that the error may not 
have been significant. In similar studies using this 
equipment the calibration precision was given as the ratio 
of the standard deviation to the mean. Using that 
criterion the results of the foregoing experiment would 
imply an uncertainty in calibration precision approaching 
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Figure 5. Worst case test for errors in calibrations. 

The envelope represents slope and intercept 
uncertainties (two standard deviations from 
mean values) . 
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8% for small loads. 

Figure 6 represents a regression line for a typical 
calibration. Note that the intercept on the ordinate is 
negative. This is the case for most of the tubes. The 
negative offset is due to the method of determining the 
area under response curves. A slight upward shift in the 
background signal from the detector occurs when the tubes 
are heated, because the warmer gas causes the acoustic 
resonance to be increased in frequency inside the detector 
cell. Therefore baselines at the end of a analytical run 
are higher than at the start. The shift in baseline 
response results in the negative offset in area. figure 1 , 
an example of a typical analysis, illustrates this 

haracteristic. 

Sampling Method 

While sampling for ammonia is relatively simple using 
this technique, precautions must be taken to ensure that 
representative samples of the ambient air are obtained. 
Contamination by nearby personnel or other sources (such as 
vehicles) is always a possibility. For this reason the 
sampling apparatus must be left largely unattended or 
situated upwind of the observer, and handling of the tubes 
must be minimized. Replicate samples help reduce the 
possibility of spurious measurements due to contamination. 
Unfortunately, for this study a shortage of collection 
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Calibration Line for Tube #6 



ammonia loading in nanograms 

Figure 6. Point plot of a typical calibration of an indi- 
vidual collection tube, accompanied by the best 
fitting line for the data. 
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tubes frequently prevented simultaneous- sampling. 
Therefore extra attention was given to avoiding 
contamination wherever possible. 

To prevent leakage of ambient air into the tubes , they 
were sealed with parafilm prior to and following each 
measurement or calibration run. Blank tubes were 
periodically checked for seepage over extended periods. 
Generally little response was noted, even after weeks of 
storage. To eliminate the possibility that samples might 
deteriorate while awaiting analysis, an experiment was 
performed to test for loss of response of a loaded tube 
after being stored overnight. The results were negative, 
indicating the tubes could be put aside for several days 
after sampling for later analysis. 

Determination of ammonia concentrations was 
accomplished by dividing the estimated load (derived for 
the relevant regression equation) by the volume of air 
pulled through the tube during the course of sampling. 
Therefore particular care was exercised in ensuring that 
the flow remained at 1 1/min through the tubes. Whenever 
it was noted that flow rates at the end of sampling varied 
by more than 5% from this value, corrections for this 
variation were included in the calculations. Sampling 
times varied a great deal, depending on the suspected 
concentration of ammonia. For extended sampling times of 
30 minutes or more periodic checks of the flow rates were 
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made. 


A flow rate of I 1/min ensured that laminar flow was 
maintained throughout the the length of the coated portion 
of the tubes (4J 130]. Laminar flow helped prevent 
particles containing aunmonium from impinging on tube walls 
and subsequently affecting the analysis. Migration of 
particles to the tube walls because of gravitational 
effects was prevented by orienting the tubes vertically 
within the mounting bracket. While particles, because of 
their inertia, were pulled directly through the tubes, 
collection efficiency for ammonia was excellent. 
Experiments made during this study confirmed that 
collection efficiency exceeded 99%, even for large loads. 

Two sites were used for sampling the ambient air. One 
of these was from the roof of the building where the 
analytical apparatus was housed (bldg. 1201 at NASA 
Langley, North latitude: 37.09 degrees and West longitude; 
76.39). Samples were taken in September 1980, January, 
February, and March 1981 were from this location, situated 
approximately eight meters above ground level. During the 
fall sampling took place two meters above the ground from a 
location one-half mile north of the original site. The 
move was made for two reasons; to avoid interfering with 
the work of other personnel using the roof, and to negate 
any effects which pollutants from a nearby steam plant 
might have had on the measurements. Transportation 
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problems necessitated the return move to the original site 
after the end of 1980. 
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VI. RESULTS AND CONCLUSIONS 
Data Redaction and Analysis 

Table I in Appendix A p; jsents the results of over 230 
measurements taken from September 1980 until mid-March 
1981. Times for initiation of sampling and estimated wind 
directions are included alongside the data. The raw data 
have been graded according to the relative accuracy of the 
given concentrations. On many days when samples were 
gathered, especially in late fall and early winter, ambient 
ammonia levels were too low to be detected by the system. 
Sampling time was extended in order to draw more air 
through the tubes and thus increase sensitivity. Yet even 
when sampling for one or two hours, frequently no ammonia 
was detected in the tubes. At other times the extended 
sampling resulted in unexpectedly large loads and caused 
off scale responses on the recorder. The former cases have 
been assigned a grade of "C." The values given represent 
one half of the minimum concentration which could be 
detected by the equipment under the circumstances. This is 
admittedly a crude estimate of the actual amount of ammonia 
present at the time, but it was felt that this information 
would be better than none at all. Those cases which are 
graded "B" correspond to off scale readings. The actual 
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loads were subsequently estimated by giving the tubes large 
calibration loads and qualitatively comparing the off scale 
responses on the recorder to those from the measurement 
runs. The estimates are better than those graded "C," 
generally representing less than 30% uncertainty in the 
accuracy of the measurements. Those measurements graded 
"A" lay within the dynamic range of the analytical 
equipment. The uncertainty intervals for these 
measurements are the sum of the uncer iiainties in 
calibrations (approximately 8% for small loads) and 
sampling flow rates (less than 5%). Therefore these values 
should differ at most from actual concentrations by 13%. 
However, a series of replicates taken in February and March 
of 1981 tend to refute this figure, since they differ form 
one another by an average of nearly 25%, suggesting that 
contamination of samples was a bigger problem than had been 
suspected. 

Table II in Appendix A gives the average daily 
concentrations when samples were taken. These have also 
been assigned grades i "A" for days when all measurements 
lay within detection limits, "B" for days when one or more 
off scale responses were noted, and ”C" for days when 
concentrations were at least once below the detection 
limits. Accompanying the measurement data are the average 
daily temperatures, wind directions, and dewpoint 
temperatures, as well as the number of samples gathered 
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each day. 

Meteorological data appearing in these tables were 
obtained from two sources. Prior to January 1981 the NASA 
Langley Energy and Environmental Monitoring System provided 
the information. Instrument failure prevented any data 
from being obtained from NASA sources after December. 
Nearby Langley Air Force Base furnished meteorological data 
during the remaintier of the study. In addition, on-site 
meteorological observations were made at sampling tiinej. 
Close inspection reveals that the furnished data differ 
little from wind and temperature estimates made at the 
measurement sites. 

Relationships with Various Meteorological Parameters 

Prom the data an obvious conclusion can be drawn; 
ammonia concentrations locally reach a minimum value in 
early winter. Figure 8 graphically depicts average daily 
concentrations during the six and one-half months of this 
study. Superimposed on this point plot are lines 
connecting the monthly means of ammonia concentrations. 
The January minimum is followed by a distinct upward trend 
in February, when some of the highest recorded 
concentrations for the entire period occurred. A possible 
explanation fnr the early winter minimum follows from an 
examination of the relationship between air temperature and 
ammonia concentrations. In table III in Appendix A daily 
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CONCENTRATIONS ppb 


Figure 8. Mean daily and monthly ammonia concentrations. 
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average concentrations have been categorized according to 
five temperature intervals and the mean value within each 
category has been presented. An upward trend in ammonia 
concentrations is discernable as temperatures increase. 
The higher mean concentration for the 50 to 59° F interval 
compared to the mean for the 60-69° F interval may be due to 
some seasonal bias in the data. Many of the days where 
temperatures averaged in the fifties were in February, a 
month of anomalously high ammonia concentrations. From the 
table there is some indication that temperatures below 50°F 
are not conducive to ammonia production. 

While table III is suggestive of some relationship 
between ammonia concentrations and air temperatures, more 
meaningful results can be obtained by means of a graphical 
presentation of the data. In figure 9 daily mean 
concentrations have been plotted against mean daily 
temperatures. A linear regression analysis has also been 
performed for this data and the resulting best fitting line 
plotted on the same graph. While ammonia concentrations 
may not be a linear function of air temperature, the 
results of the analysis point to some relationship. In 
fact a linear correlation coefficient of 0.5346 for these 
data indicates that there is a 99.9% probability that the 
dependent and independent variables are correlated [3J. In 
light of these findings a partial explanation for the 
upward trend in February concentrations may be inferred if 
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the meteorological record for that month is examined 
closely. February 190i was notable for several abnormally 
warm periods following unusually cold weather in January 
and late December. The mean temperature for February was 
43.1°F (6.2°C), compared to 32.6°F (0.3°C) in January and 
41. 6° F (5.3°C) in December. 

Scientists are acutely aware of the danger of 
inferring a cause-and-ef feet relationship between two 
parameters simply because they appear correlated. If 
ammonia levels were a function of air temperature alone, 
then measured ammonia amounts would have been significantly 
higher in October and November 1980 (mean temperatures: 
61.5° F and 49.1 °F, respectively) than in February. Other 
factors must be considered as well as air temperature. 

Hoell et a_l. [20] noted a similar increase in local 
ammonia concentrations in March of 19 79, when levels a.^, 
high as 10 ppb were measured. It was tentatively concluded 
that volatilization of gaseous ammonia from nearby fields 
where fertilizers had been recently applied largely 
contributed to the March maximum in that year. An attempt 
was made for this study to further investigate the 
prevailing agricultural practices of surrounding stt-as 
regarding fertilizer application. Several calls 
placed to agricultural agents in nearby Virginia 
localities. Mr. Jim Belote [2], an agricultural extension 
agent for Virginia Beach, stated that most fertilizers were 
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applied in April through May. However# a top dressing of 
fertilizer (approximately 80 lbs. per acre) was applied to 
winter wheat beginning in late Peb)r«sr/. He estimated that 
13#000 acres are given to winter wheat production in 
Virginia Beach. Mr. Ben S. Lee [26] in Southhampton 
County stated that some farmers applied fertilizer in fall, 
while most waited until March. But he also indicated that 
large amounts of fertilizer wer- applied to winter wheat 
crops in February. Mr. Lee estimated that the total 
acreage in winter wheat was nearly 10,000 acres in 
Southhampton County. While other extension agents could 
not be reached, office personnel in several county 
agricultural extension offices agreed that significant 
amounts of fertilizer are generally applied to winter wheat 
crops, weather permitting, beginning in February. 

Since most of the fertilization of winter wheat in 
surrounding areas appears to have taken place after high 
ammonia concentrations first appeared in early February, 
positive conclusions cannot be drawn about about any link 
between the two events. It may be useful to examine other 
factors besides temperature or agricultural practices that 
could have affected local ammonia production in February 
1981. 


Dawson 

[8] 

cited soil 

moisture as 

well as soil 

temperature 

as 

an 

important 

parameter 

regulating the 

release of 

ammonia 

from the 

soil by 

microb iolog ical 
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activity. He indicated that saturated soils release little 
ammonia, while completely dry soils cannot support the 
microorganisms which convert organic substances into 
gaseous ammonia. Intermediate levels of soil moisture are 
required for optimal volatilization of ammonia from the 
soil. The best method for determination of the water 
content of a particular soil is to continuously monitor 
moisture levels. Unfortunately such data was not available 
for this research. Instead a water budget model which 
estimates soil moisture based on evapotranspiration rates 
has been used to determine whether soil moisture levels may 
have varied significantly during the course of this study. 
The mo^lel requires information on soil field capacity (a 
measure of how much water the soil can hold) , monthly mean 
temperatures, and monthly rainfall amounts in order to 
arrive at estimates of soil water content. When soil field 
capacity is reached (80% of saturation, according to 
Dawson) any excess rainfall is categorized as runoff. 
Field capacity for this case was estimated by referring to 
soil surveys for nearby counties [43] . Research has shown 
that the model successfully predicts actual water runoff 
[24]. However, no attempt was made to verify the model's 
estimates of soil moisture. Results of the model siiould 
therefore not be interpreted too narrowly. 

For the fall to late winter period covering this 
study, the model showed that water levels rose from near 
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zero in late September to 25% of field capacity by November 
first. In contrast moisture was nearly 30% of field 
capacity by the end of February. Dawson estimated that 
ammonia volatilization peaks at moisture content 
corresponding to 20% of field capacity. Model results tend 
to refute the hypothesis that the February ammonia levels 
were due in part to soil moisture. 

The water budget model may not be very useful for 
interpreting February's relatively high ammonia 
concentrations. But it can be applied with more conclusive 
results when comparing the data presented here for 
September 1980 with similar measurements for September 1979 
published by Hoell ^ al . [20]. Ammonia levels in September 
1979 averaged around 2 ppb, over twice as high as those 
found in September 1980 at the same location. The model 
shows that extensive rainfall throughout the summer of 1979 
resulted in soil moisture remaining at nearly 57% of field 
capacity for September cf that year. By contrast severe 
drought in 1980 completely depleted the soil of available 
moisture by July. Dry conditions continued until 
moderating temperatures and more plentiful rainfall brought 
soil moisture levels back up in the fall. The drought must 
have had a significant effect on microbiological activity 
in local soils. (See Appendix B for model output.) 

Results presented here confirm the importance of air 
temperature and rainfall (and therefore soil temperature 
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and moiiiture) in regulating ammonia production and ambient 
ammonia concentations. While rainfall tends to increase 
ammonia production from previously dry soils, it also may 
scavenge large amounts of ammonia from the air. Georgii 
and Muller [11] witnessed this effect, as did Harward £t aL 
[14]. Results from this study, however, proved 
inconclusive. On several occasions ammonia concentrations 
indeed dropped following rain; on others concentrations 
apparently increased. Insufficient data also inhibit the 
formation of any conclusions. Lack of a suitable shelter 
prevented any sampling during rainfall. Furthermore, since 
most rainstorms during the study were not showers, but 
prolonged rainfalls, follow-up sampling during normal work 
hours was frequently prevented. 

Other meteorological parameters besides temperature or 
rainfall may have an impact on ambient ammonia 
concentrations. In a theoretical study Lau and Charlson 
[25] concluded that high humidities may reduce the amount 
of gaseous ammonia in the air by promoting its conversion 
to dissolved ammonium sulfate in haze or cloud droplets. 
Data on dewpoint temperatures were made available for this 
study. Therefore the average water vapor mixing ratios and 
relative humidities could be calculated, and both were 
correlated with mean ammonia concentrations as was done for 
air temperature. The results appear in graphical form in 
figures 10 and 11. 
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Figure 11. Mean daily concentrations plotted against mean 
daily relative humidities, accompanied by the 
best fitting line for the data. 
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Contrary to the assertion of Lau and Charlson, ammonia 
levels seemed to actually increase with higher relative 
huftiidities. The correlation coeficient for the regression 
line in figure 11 is approximately 0.367, which for N«47 
implies a high probability (greater than 95%) that relative 
humidity and ammonia concentration were directly 
correlated. Even more interesting is the high correlation 
between the water vapor mixing ratio of the air and ambient 
ammonia levels. The linear correlation coefficient here 
(0.722) was much higher than those for either relative 
humidity versus concentration or temperature versus 
concentration. This relationship was not due to 
measurement errors caused by collection of water on the 
sampling tubes. Experiments were performed on the tubes' 
sensitivity to various humidities, both during this study 
and others [14] , and results were negative. 

Although explanations for the strong relationship 
between water vapor mixing ratio and ammonia concentration 
must remain in the realm of speculation, an hypothesis can 
be advanced here. Ammonia levels seem to be related to air 
temperature and relative humidity separately. When 
correlating water vapor mixing ratio and ammonia 
concentrations, we are essentially adding these two 
effects. In other words, if ambient ammonia tends to 
increase as either the relative humidity or air temperature 
increases, then it should certainly be high under warm. 
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humid conditions. This still does not explain the 
relationship betvt’een relative humidity and ambient ammonia 
levels. One clue may come from the frequently observed 
high variations in ammonia content from one sample to the 
next. Ammonia often seems to arrive in "puffs" or events, 
rather than remaining at similar concentrations. On humid 
days when clouds are present, evaporation at cloud bases 
creates "dry" aerosols from dissolved species (consisting 
of at least some ammonium nitrate) in the cloud droplets. 
Research has shown that the equilibrium chemistry of 
ammonium nitrate is highly temperature dependent. Warming 
of the ambient air (for exaiiiple, by adiabatic descent) 
results in formation of gaseous ammonia and nitric acid 
from the dissociation of the solid ammonium nitrate [4]. 
Thus "clouds" of gaseous ammonia and nitric acid may 
frequently descend to the surface from the condensation 
level. This effect would be most prevalent on warm days 
when vigorous mixing takes place. Indeed, some of the 
highest concentrations were noted on partly cloudy days. 

No quantitative data can be presented here to further 
examine the high variations in ammonia content within an 
individual air mass. Such a study would require an 
extensive monitoring network. But from a single station 
different air masses can at least be crudely categorized 
according to ammonia levels. This is usually accomplished 
by means of concentration wind roses. The sampling 
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locations for this research offer an opportunity to 
investigate the contention that most atmospheric ammonia is 
of continental origin. If the continents constitute the 
major source for atmospheric ammonia, air masses having 
long trajectories over open water would therefore contain 
less ammonia than those which generally passed over land 
before reaching coastal Southeastern Virginia. Table IV 
(in Appendix A) presents the mean ammonia concentrations 
for individual measurements classified according to 
prevailing wind directions at the time of sampling. Figure 
12 illustrates these results by means of a concentration 
wind rose superimposed on a map of the area. Ammonia 
concentrations were significantly higher when winds were 
from the south or southwest than when winds were from the 
east, northeast, or north. However, air masses arriving 
from the west or northwest did not contain significantly 
more gaseous ammonia than those coming off open water. 
Furthermore, highest concentrations were found when 
southeasterly winds prevailed. This anomalous maximum was 
probably due to anthropogenic sources. In particular, an 
industrial source for gaseous ammonia — a fertilizer 
plant — is situated roughly 10 km southeast of NASA Langley. 
Air arriving from the other seven wind directions was 
comparatively pristine. Of course surronding areas are by 
no means undeveloped, and the least development has 
occurred to the north and west of Langley. This may 
partially explain the comparatively low ammonia content 


Figure 12. Concentration wind rose for individual 
measurements . 
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found in the air arriving from the west, northwest, or 
north. The humidity factors mentioned above may also have 
influenced the results. Cool temperatures and low 
humidities generally prevailed under north or northwesterly 
winds, while winds from the south or southwest were usually 
accompanied by warmer and moister conditions. A further 
biasing of the data may be inferred when examining those 
days when easterly or northeasterly winds prevailed. Most 
of these cases occurred in the early fall, when ammonia 
levels were generally higher overall. By contrast the 
winter sampling was characterised by an overabundance of 
northerly and northwesterly winds. 

Because of these reservations about the available data 
base, any conclusions about the sources for most of the 
"natural" ammonia in the air of Southeastern Virginia must 
remain tentative. Yet continental air masses arriving from 
the south and southwest did contain significantly more 
gaseous ammonia than maritime air masses. Man's 
participation in the production of atmospheric ammonia, at 
least around this locality, was also quite apparent from 
the data. 

Implications of Results 

One of the major problems encountered during this 
research involved selecting appropriate sampling intervals 
and adjusting the sensitivity of the analytical apparatus 
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so that ammonia could be measured accurately. Large, 
short-term variations in ammonia concentrations were 
frequently encountered that exceeded the dynamic range of 
the equipment. This demonstrates that gaseous ammonia 
simply cannot be thought of or treated as a uniformly 
distributed background constituent. This is especially 
true near the ground, the source of virtually all gaseous 
ammonia. When numerical modelers talk of steady-state 
conditions in the atmosphere, they almost always are 
oversimplifying actual conditions. Of course they are 
acutely aware of this, and would respond that there is 
usually no better way to analyze atmospheric behavior, 
whether it be of a dynamical or chemical nature. Ammonia's 
heterogenous chemical interactions with other atmospheric 
constituents present an almost intractable problem even 
when steady-state conditions are assumed. Once these 
problems are solved, scientists can move on to more 
realistic assumptions. At that time they will require 
accurate information about the real-time variations in 
concentrations of a host of trace species in the air. For 
measurement of ammonia, at least, such capabilities simply 
do not exist. Integral results such as those from tliis 
technique may be adequate for now, but a reliable 
fast-response monitoring system with a large dynamic range 
in sensitivity must be developed. 
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Besides finding that concentrations of ammonia vary a 
great deal, th^s study has shown that ammonia levels in 
southeastern Virginia are much lower than those reported 
from other measurement programs in both Europe and other 
places on this continent. Some of the higher measurements 
given for other areas may reflect inaccurate measurement 
techniques. For example, some have cited evidence that 
measurement errors with wet chemical techniques may 
approach 1 ppb (25]. On the other hand, the measurements 
generally found in this study do have a physical 
explanation (which hinges on soil conditions). 

Junge (21], when examining rainwater analyses made in 
this country, noticed that dissolved ammonium in rainwater 
was much less (sometimes almost absent) across the entire 
Southeast than elsewhere in the United St£tes. He 
attributed this to systematically lower gaseous ammonia 
concentrations in the air, and hypothesized that some 
factor may prevent the escape of ammonia from soils in the 
Southeastern U.S.A. He found that a yellow-red lateritic 
soil type is predominant throughout the Southeast. Soils 
in the area have an average pH of less than six (21]. 
Soils of low pH inhibit both microbiological activity and 
the release of gaseous ammonia into the air. Coastal 
Virginia also has markedly acidic soils. According to soil 
surveys (43] pH values of less than 4.5 prevail for most 
soils in the Dismal Swamp region, while soils in much of 
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the rest of Tidewater Virginia have pH values varying from 
3.6 to 6< 0 • 

Since natural emissions of gaseous ammonia may be 
inhibited by soil acidity/ anthropogenic sources are likely 
to be more noticeable locally than where background levels 
of ammonia are higher. Such was the case with this study. 
A small fertilizer plant situated miles from the sampling 
locations apparently produced some of the highest 
concentrations that were encountered. Systemat ica*..y 
higher levels of ammonia in late winter may have been 
caused by man's intervention in natural soil processes. By 
^PPlyirig fertilizers to fields farmers directly introduce 
ammonia into the air by increasing volatilization of 
ammonia from soils saturated with nitrogen. They may 
indirectly affect long-term emissions of ammonia by 
altering the pH of the soil through fertilization. 
Efficient and economical use of limited resources, as well 
as concern for possible effects of this vast dumping of 
nitrogen compounds into the air, should lead 
agriculturalists to closely examine existing practices with 
an eye toward reducing ammonia emissions from agricultural 
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Comments on the Efficacy of This Technique for Routine 
Monitoring 

It was unfortunate that more could not be iiaid in this 
paper about short-term variations of ammonia content in the 
atmosphere. It was hoped that links could be established 
between ammonia concentrations and specific meteorological 
events, such as frontal passages or rainstorms. Such was 
not the case, however, and some of the shortcoming in the 
data must be attributed to the measurement technique 
itself. Although LPS is far superior to most other 
techniques of measurement, some of its limitations and 
drawbacks should be commented on briefly. 

It is apparent that low ammonia concentrations (as low 
as a few tens of parts per trillion by volume at times) 
make real-time measurements with existing technology 
extremely difficult. Integral methods employing 
preconcentration of some kind are the only available 
alternative. Yet some refinements to the system are quite 
possible without significantly altering the basic 
technique. Automated sampling would be an important 
improvement. As the disappointing results of replicate 
sampling have shown, contamination by contact with the 
collection tubes may constitute a major problem. Poor 
replicate results may also stem from inconsistencies in the 
collection tubes' performance. Early in .this study a 
gradual loss of sensitivity for some of the tubes was 
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noticed. This necessitated frequent and time-consuming 
recalibration of the tubes. Eventually some tubes' lack of 
sensitivity or reproducibility made re juventation by 
oxidation necessary. This wasted further time and did not 
always work. Eventually serendipity led to the discovery 
that a simple washing with distilled water completely 
restored a tube's performance. This was especially 
fortunate because longer sampling times had made it 
necessary to rejuvenate the tubes more often. 

More distressing than gradual aging of the tubes was 
their tendency to occasionally become "poisoned" by the 
action of some unknown agent. Rejuvenation restored them 
but did not restore confidence in the accuracy of 
measurements taken with poisoned tubes. This may explain 
some of the poor reproducibility noted in replicate 
sampling. The possibility of poisoning also made constant 
monitoring of each tube's performance an absolute 
necessity, which wasted still mors precious time better 
spent samplin';^. The exact cause for sudden loss of 
sensitivity has not been determined. It is suspected, 
however, that sulfur compounds may have contributed to this 
phenomenon. A steam generating plant was located only a 
few hundred meters from the original sampling site. Sudden 
poisoning was more frequently encountered when winds 
carried the plume containing large amounts of sulfur 
dioxtide from this plant toward the sampling station. For 
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this reason the tubes' usefulness in extremely polluted 
environments must be questioned. 

It is possible to envision the entire LPS system being 
incorporated into a single monitoring device no larger than 
conventional monitors for other trace gases. The 
technology is presently available, but modifications such 
as a smaller laser and the inclusion of a computer to allow 
automated analysis would require considerable ingenuity and 
technical expertise, not to mention large sums of money. 
The overall performance of the system was encouraging 
enough to merit its further development. There remains the 
possibility of incorporating basic design changes into the 
system to permit real-time measurements of ammonia. 
Perhaps interferometry could be used to detect minute 
changes in the index of refraction within the gas cell as 
ammonia molecules absorb radiation and heat the medium. 
Eventually an attempt must be made to bypass the currently 
necessary step of preconcentration. Fast response in 
real-time monitoring is the goal of any ^ situ measurement 
technique . 
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Tables of Data 
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TABLE I. 

INDIVIDUAL AMMONIA MEASUREMENTS 

Date Time Concentration Wind Grade 

(ppb) Direction 


09-03-80 

9:35 

0.20 

N 

C 


9:58 

0.20 

N 

G 


11:55 

0.20 

N 

C 


12:56 

0.66 

NE 

A 


1:20 

0.20 

NE 

C 


1:42 

0.20 

NE 

C 


2:05 

1.08 

E 

A 

09-0^-80 

10:26 

1.68 

SE 

A 


11:48 

2.22 

SE 

A 


12 ;33 

1.26 

SE 

A 


1:17 

1. 18 

SE 

A 


1:58 

1.98 

SE 

A 

09-08-80 

9:50 

1.50 

sw 

A 


10 :22 

1.32 

NW 

A 


10 :53 

1.16 

NW 

A 


12:22 

0.20 

N 

C 


12:54 

1.10 

NE 

A 


1:26 

1.08 

NE 

A 


1:57 

1.56 

NE 

A 


2:29 

0.20 

NE 

G 


3:00 

0.20 

NE 

C 

09-10-80 

8:40 

0.80 

SW 

A 


9:49 

0.88 

SW 

A 


10:29 

0.20 

SW 

C 

09-19-80 

9:26 

1.26 

N 

A 


10:11 

1.12 

N 

A 


12:23 

0.96 

N 

A 


1:09 

0.86 

NE 

A 


ls56 

1.12 

NE 

A 


2:40 

0.80 

NE 

A 

09-22-80 

9:30 

0.90 

SW 

A 


9:52 

2.44 

SW 

A 


10:20 

0.90 

SW 

A 


10:12 

1.20 

SW 

A 


11:11 

1.20 

SW 

A 


11:32 

1.62 

SW 

A 


11:54 

1.08 

SW 

A 


1:08 

1.02 

VJ 

A 


1:49 

0.62 

W 

A 


2:28 

1.08 

w 

A 


3:12 

1.30 

w 

A 

09-^4-90 

9 : 05 

0772 

N 

A 


9:51 

0.46 

N 

A 


10 : 36 

0.62 

NE 
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TABLE I. 

CONTINUED 

Date Time Concentration wind Grade 

(ppb) Direction 





NE 
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12:34 

0.76 

NE 

A 


1:15 

0.48 

NE 

A 


1:45 

0.62 

NE 
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2:31 

0.70 

NE 

A 


3:13 

0.66 

NE 

A 

09-i6-80 



SW 

A 


9:20 

1.08 

SW 

A 


10:00 

0.86 

W 

A 


10:40 

0.80 

NW 
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12:20 

1. 10 

NW 

A 


1:00 

0.94 

NW 

A 


1:40 

0.90 

W 

A 


2:20 

0.62 

NW 

A 


3:00 

1.40 

N 

A 

fiI;gfctgf:TOMi 

$;20 

5TT? 

NE 
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0.18 

E 

A 


10:40 

0.40 

E 

A 


12:30 

0.36 

NE 

A 


2:00 

1.02 

E 
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2:50 

0.70 

E 

A 


3:30 

0.10 

E 

A 

iVtZSBEEFIiiHH 


oi 

NW 

C 


12:45 

0. 10 

NW 

C 


1:10 

0.36 

mi 

A 


1:35 

0.27 

NW 

A 

1o-08«§5 

lO : 4S 

oTIB 

s 

A 


11:10 

0.55 

s 

A 


12:17 

0.48 

s 

A 


12:40 

0.56 

s 

A 

10-10-80 

TTT5 

0TT5 

NE 
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9:02 

0.05 

NE 
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9:25 

0.10 

NE 
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9:47 

0. 10 

NE 
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10:09 

0.04 

NE 

c 

10-15-80 
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9:52 

0.05 
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1:16 

0.82 

s 
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1:39 

1.03 
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2:01 

1. 32 
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oTo3 

NW 
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1:51 
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2:19 
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C 

10-22-80 

9:50 

0.50 
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TABLE I. 
CONTINUED 


Date Time Concentration Wind Grade 

(ppb) Direction 
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TABLE I. 
CONTINUED 


Date 


Time Concentration Wind 

(ppb) Direction 


Grade 
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TABLE I. 
CONTINUED 


Date 


Time 


Concentration 

(ppb) 


Wind 

Direction 


Grade 



TABLE II 


MEAN DAILY CONCENTRATIONS AND 



Ave. 

Ave. 

Ave. 

Date 

Concen. 

Temp. 

Wind 


(ppb) 

( P) 

Dir. 

09-03-80 

0.48 

80 

NE 

09-05-80 

1.66 

78 

SE 

09-08 -80 

1.00 

74 

NE 

09-10-80 

0.70 

72 

SW 

09-19-80 

X.02 

74 

N 

09-22-80 

1.21 

84 

SW 

09-24-80 

0.64 

73 

NE 

09-26-80 

0.99 

75 

W 

09-29-80 

0.42 

68 

E 

10-06-80 

0.19 

59 


10-08-80 

0.51 

65 

S 

lO-lO-SO 

0 . 10 

69 

E 

10-15-80 

0.6S 

61 

S 

10-20-80 

0.06 

58 

NW 

10-22-80 

0.17 

62 

NW 

10-29-80 

0.08 

50 

N 

10-31-80 

0.02 

50 

W 

11-03-80 

0.11 

50 

E 

11-07-80 

0.20 

57 

SW 

11-10-80 

0.60 

59 

W 

11-14-80 

0.09 

60 

SW 

11-19-80 

0.09 

40 

NW 

11-21-80 

0.09 

43 

N 

11-24-80 

1.00 

58 

SE 

11-26-80 

0.05 

41 

N 

li-01-&0 

oTI^ 

56 

SW 

12-05-80 

0.06 

41 

NW 

12-08-80 

0.26 

63 

SW 

12-12-80 

0.12 

44 

w 

12-15-80 

0.04 

43 

E 

12-17-80 

0.05 

34 

N 

ol-li-8i 

■jjDHjfefljjl 

T5 ” 

NW 

01-14-81 

0.07 

32 

SW 

01-15-81 

0.20 

36 

W 

01-16-81 

0.06 

34 

SE 

01-19-81 

0.11 

40 

S 

01-23-81 

0.05 

43 

W 

02-0^-lT“ 

TToo 

50 

W 

02-03-81 

0.29 

26 

w 

02-04-f 1 

0.22 

25 

w 

02-09-81 

0.15 

36 

w 

02-11-81 

2.21 

56 

SE 

02-13-81 

0.01 

30 

NE 


“■"4 
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1 


1 


METEOROLOGICAL DATA 


Ave. 

No. 


Dewpoint 

of 

Grade 

( F) 

Runs 



60 

7 

C 

61 

5 

A 

49 

9 

C 

53 

3 

C 

55 

6 

A 

61 

11 

A 

52 

9 

A 

49 

9 

A 

44 

7 

C 

32 

4 

C 

39 

4 

A 

52 

5 

C 

34 

5 

C 

33 

3 

C 

37 

6 

C 

26 

2 

A 

2 7 

4 

C 

27 

4 

G 

28 

4 

A 

31 

3 

A 

26 

4 

C 

13 

2 

C 

21 

2 

C 

39 

1 

B 

14 

3 

C 

22 

3 

C 

10 

3 

C 

33 

3 

A 

17 

1 

A 

19 

3 

C 

15 

2 

c 

2 

2 

A 

18 

2 

C 

30 

1 

A 

31 

1 

A 

22 

4 

A 

31 

2 

A 

39 

5 

C 

7 

4 

C 

8 

5 

C 

16 

4 

A 

52 

11 

A 

16 

5 

C 
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TABLE II. 
CONTINUED 


Date 

Ave. 

Concen. 

(ppb) 

Ave. 
Temp. 
( F) 

Ave. 

Wind 

Dir. 

Ave. 

Dewpoint 
( F) 

No. 

Of 

Runs 

Grade 

S2- 18-81 

1. 36 

56 

S 

48 

8 

A 

02-20-81 

0.37 

53 

W 

50 

8 

A 

02-23-81 

0.38 

55 

S 

46 

8 

G 

02-25-81 

0.81 

45 

SW 

45 

10 

A 

05-69-61 

5712 

NA 

NA 

NA 

6 

A 

03-13-81 

0.93 

NA 

NA 

NA 

2 

A 
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TABLE III. 

MEAN DAILY CONCENTRATIONS ACCORDING TO TEMPERATURE 

INTERVALS 


Mean 

Mean 

Number 

Temperature {°F) 

Concentration (ppb) 

of Days 


< 40 

0.12 

9 

40 - 49 

0.16 

9 

50 - 59 

0.56 

14 

60 - 69 

0.31 

7 

70 

0.96 

8 
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TABLE IV. 

CONCENTRATION WIND ROSE BY INDIVIDUAL OBSERVATIONS 

Number 

Wind Direction Mean Concentrations of Observations 


N 0.30 26 

NE 0.38 25 

E 0.31 13 

SE 1.61 20 

S 0.81 22 

SW 0.63 51 

W 0.40 46 

NW 0.37 20 
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APPENDIX B 

Program THORN and Its Output 

Appendix B C’^'^ntains the results of the water budget 
model used in this study to estimate the water content of 
the soil. For purposes of comparison the model has been 
run for a five year period from 1977 to 1981. Data for the 
months after February of 1981 were of course not available, 
but dummy data were inserted to allow the model to run 
through the entire year. Folowing the output is a listing 
of the FORTRAN code comprising the model. 
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TDEGPaTEMPERATURE IN DEGREES P 
TDEGC-TEMPERATURE IN DEGREES C 
PPTIN-MONTHLY PRECIPITATION IN INCHES 
PPTCM-MONTHLY PRECIPITATION IN CENTIMETERS 
SMALLI»MONTHLY HEAT INDEX 

UPETCM-UHADJUSTED POTENTIAL EVAPOTRANSPIRATION IN CM 
CORPHI-CORRECTION PACTOR POR LENGTH OP DAYLIGHT AT STATION 
APETCMa ADJUSTED PET (»CORPHI*UPETCM) 
RMNDRaPRECIPITATION-ADJUSTED PET 
STRCHGaSTORAGE CHANGE 
STRAGEaWATER IN SOIL (CM) 

%CAP»PER CENT OP SOIL PIELD CAPACITY 

WTRDEPaWATER DEPICIT (NO MORE WATER IS STORED IN SOIL) 
WTRSPLaSURPLUS OP WATER BEYOND PIELD CAPACITY 
ACTEVP* ACTUAL EVAPOTRANSPIRATION IN CM 
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C THIS PROGRAM COMPUTES POTENTIAL 
EVAPOTRANSPIRATION (USING 

C THE THORNTHWAITE METHOD), THEN COMPUTES 
ACTUAL EVAPOTRANSPIRATION 
C MONTHLY AT A GIVEN SITE. 

C ♦ 

C PLEASE NOTE; AT LEAST TWO DATA FILES ARE 
REQUIRED TO RUN THIS 

C PROGRAM; DAYDEC.DAT AND A DATA FILE 
CONTAINING MONTHLY TEMPERATURES 
C AND PRECIPITATION AMOUNTS. 

C 

DIMENSION AVMDEC(24), DAY.SM(24), 
TDEGF(24), TDEGC{24), 

1 PPTIN(24), PPTCM(24), SMAI.LI(24), 
UPETCM{24) ,CORPHI(24) , 

2 APETCM(24), RMNDR(24), STRCHG(24), 
STRAGE(24), WTRDEF(24), 

3 WTRSPL(24), SNAME(8), ACTEVP(24) 
DIMENSION FNAMEd ) ,PCMAX(24) 

DOUBLE PRECISION FNAME 

C AVMDEC-AVERAGE MONTHLY SOLAR DECLINATION 
FOR THE STATION 

C DAYSM»NUMBER OF DAYS IN A MONTH 

C TDEGF-TEMPERATURE IN DEGREES FARHENHEIT 

C TDEGC-TEMPERATURE IN DEGREES CENTIGRADE 

C PPTIN=MONTHLY PRECIPITATION IN INCHES 

C PPTCM=MONTHLY PRECIPITATION IN 

CENTIMETERS 

C SMALLI=MONTHLY HEAT INDEX 

C UPETCM=UNADJUSTED POTENTIAL 

EVAPOTRANSPIRATION IN CENTIMETERS 
C CORPHI?»CORRECTION FACTOR FOR LENGTH OF 

DAYLIGHT AT STATION 

C APETCM=ADJUSTED PET (sCORPHI*UPETCM) 

C RMNDR=PRECIPITATION- ADJUSTED PET 

C STRAGEs»WATER IN THE SOIL (CM) 

C STRCHG=STORAGE CHANGE 

C WTRDEF«WATER DEFICIT (NO MORE WATER IS 

STORED IN SOIL). 

C WTRSPL=SURPLUS OF WATER BEYOND TCE 

MAXIMUM STORAGE CAPACITY 

C ACTEVP=THE ACTUAL EVAPOTRANSPIRATION IN 

CM 

OPEN(UNIT=21 ,ACCESS='SEQIN’ ,FILE='DAYDEC.DAT' ) 
READ (21,300) ( AVMDEC ( I ) , 1=1 , 1 2 ) 

READ (21,400) ( DAYSM( I ) , 1=1 , 1 2) 

300 FORMAT (12F6.2) 

400 FORMAT (12F5.1) 

CLOSE (UNIT=21) 

962 WRITE(5,789) 

789 FORMAT (///,' TYPE NAME OP DATA FILE 


(IF YOU WANT TO STOP 
ITYPE NONE) : ' ) 

READ (5,796) PNAME(l) 

796 FORMAT (AlO) 

IF (FNAMEd ) .EQ. 'NONE* )RETURN 

OPEN(UNIT»21 ,ACCESS»'SEQIN' ,FILE»FNAME( 1 ) ) 
PI-3.14159 
PIRAD-PI/180. 

DO 200 1-1,12 
AVMDEC( I)»AVMDEC( I) *PIRAD 
C AVERAGE MONTHLY DECLINATION IS NOW IN 
RADIANS! 

200 CONTINUE 

P-l./PIRAD 
READ (21,936) lYEARS 
READ (21,101) NY 
101 FORMAT (13) 

936 FORMAT (14) 

DO 4 1-1,12 
TDEGF(I)=0. 

TDEGC(I)=0. 

PPTIN(1)»0. 

PPTCM( I)=0. 

SMALLI(I)»0. 

UPETCM(I)=0. 

APETCM(I)=0. 

CORPHI(I)=0.0 

RMNDR(I)=0. 

STRCHG(I)-0. 

STRAGE(I)»0>0 

WTRDEP(I)=0. 

WTRSPL(I)=0. 

ACTEVP(I)=0. 

4 CONTINUE 
SUMI-0. 

SUMTF-0 . 

SUMPIN-0. 

SUMTC-0. 

SUMPCM-0. 

SUMAPE-0. 

SUMDEF-0. 

SUMSPL-0. 

SUMACT-0. 

READ (21,2) SNAME, SLAT, NORS , SLONG, 
EORW, SELEV, Z, WHC, UNIT 
2 

PORMAT(8A4,F7. 2,A2,F7. 2,A2,F7. 2 , A2 , F 7. 2 ,A2 ) 
STRAGE( 2)=WHC 
DO 102 JNY=1,NY 
IF (JNY.GT.l) GO TO 601 
WRITE (5,99) 

SNAME , S LAT , NORS , S LONG , EORW , SELEV , Z , WHC , UN I T , I Y 


EARS 

99 FORMAT (///, IX, 8A4 ,//, IX LATITUDE ' , 
F7. 2, A2,/, IX, 'LONGITUDES 
1 

F 7. 2, A2, /, 1 X, • ELEVATION • ,F7. 2, A2,/, IX, 'WATER 
HOLDING 

2 CAPACITY', F7.2, A2,/,' YEAR 
M4,///) 

GO TO 602 

601 WRITE(5,603)IYEARS 

603 FORMAT(///,36X,I4) 

602 JI-I 

IF (JNY.EQ.l) JI-3 

C THIS SOLVES AN INITIALIZATION PROBLEM 
WITHIN THE PROGRAM. 

C IN THE SECOND MONTH STORAGE SHOULD BE AT 
CAPACITY. 

READ (21,5) (TDEGF(I), I»l,12) 

READ (21,753) ( PPTIN ( I ) , I»1 , 1 2 ) 

5 FORMAT (12F 5.2) 

753 FORMAT (12P 6.3) 

DO 999 1-1,12 
PPTCM(I)»PPTIN(I) *2.54 
TDEGC( I)-( 5./9. ) *(TDEGF( I)-32. ) 

999 CONTINUE 

SUMI-0. 0 
DO 8 I-JI,12 

IF (TDEGC(I) .LE.O.) GO TO 6 
SMALLI(I)=(TDEGC(I)/5. ) **1.51 4 
GO TO 7 

6 SMALLI(I)»0. 

7 SUMI-SUMI+SMALLI (I) 

8 CONTINUE 

A»(0. 000000675* (SUMI**3. ) )-( 0. 0000771 *( SUMI* *2 

.)) 

1 +(0.01792*SUMI)+0. 49239 

C THE ABOVE CALCULATES THE STATION 

CONSTANT FROM THE 

C SUM OP THE HEAT CONSTANTS OF EACH 

STATION 

DO 10 I-JI,12 

IF (TDEGC( I ) .LE. 0. ) GO TO 9 

UPETCM( I)*l . 6*( (10.*TDEGC(I)/SUMI)**A) 

C THE ABOVE CALCULATES THE UNADJUSTED 

POTENTIAL EVAPOTRANS PI RATION 
GO TO 10 

9 UPETCM(I)=0. 

10 CONTINUE 

PHIRAD=SLAT*PIRAD 

C 

C DO LOOP 11 CALCULATES THE CORRECTION 


fact(;r from latitude, 

C LENGTH OF DAY, AND SOLAR DECLINATION 

AND THEN 

C CALCULATES THE ADJUSTED VALUES FOR 

POT EVAPOTRANSPIRATION 
C 

DO 11 I-JI,12 
Q-DAYSM(I) 

R-AVMDEC( I) 

S— S I N ( PH I P-AD ) /COS ( PH I RAD ) 
T»SIN(R)/COS{R} 

U«S*T 

V»ACOS(U) 

CORPHI(I)-(Q*V^P)/2700. 

APETCM ( I ) -UPETCM ( I ) *CORPHI ( I ) 

1 1 CONTINUE 

IF (JNY.GT.l)GO TO 80 
C 

C DO LOOP 18 CALCULATES THE MONTHLY 

WATER BUDGET FOR THE 
C STATIONS FOR THE FIRST YEAR 

C 

DO 18 I»JI,12 
STRMAX»«WHC 

RMNDR ( I ) »PPTCM ( I ) -APETCM { I ) 

12 IF (RMNDR(I).GE.O. )GO TO 1 3 

IF (STRAGE(I-1 ) .GT.O. )GO TO 13 
STRCHG(I)=0. 

GOTO 14 

13 STRCHG(I)»RMNDR(I) 

14 STRAGE(I)»STRAGE(I-1 )+STRCHG(I) 

IF (STRAGE(I) .GT.STRMAX) GOTO 15 
WTRSPL(I)«0. 

GOTO 16 

15 WTRSPL(I)=STRAGE(I)-STRMAX 
STRAGE { I ) »WHC 

STRCHG( I)«STRCHG( I)-WTRSPL(I) 

16 IP (STRAGE(I) .LB.O. ) GOTO 17 
WTRDEF(I)=0. 

GOTO 18 

17 WTRDEF(I)=RMNDR(I)+STRAGE(I-1 ) 
STRAGE(I)=0. 

STRCHG ( I ) =STRCHG ( I ) -WTRDEF ( I ) 

IF 

(STRAGE(I-1 ) .GT.O. . AND. RMNDR( I ) . GT. 0. ) GOTO 
18 

STRCHG(I)=-STRAGE(I-1 ) 

WTRDEF ( I ) =RMNDR ( I ) -STRCHG ( I ) 

1 8 CONTINUE 
GOTO 118 

80 DO 118 1=1,12 

STRAGE (0)=0.0 

RMNDR ( I ) =PPTCM ( 1 ) -APETCM ( I ) 


112 IF (RMNDRUNGE.O. ) GOTO 113 

IP (STRAGE(I~1 ) .GT.O. ) GOTO 113 
STRCHG(I)-0. 

GOTO 114 

113 STRCHG(I)»RMNDR(I) 

114 STRAGE(I)-STRAGE(I-1 )+STRCHG(I) 

IF (I.GT.l) GOTO 224 

STRAGE ( I ) -STRAGE (12) -fSTRCHG ( 1 ) 

224 IP (STRAGE(I) .GT.STRMAX) GOTO 115 

WTRSPL(I)-0.0 
GOTO 116 

115 WTRSPL(I)-STRAGE(I)-STRMAX 
STRAGE ( I )-WHC 

STRCHG ( I ) -STRCHG ( I ) -WTRSPL ( I ) 

116 IF (STRAGE(I) .LE. 0. ) GOTO 117 
WTRDEF(I)-0. 

GOTO 118 

117 WTRDEF( I)»RMNDR(I)+STRAGE(I-1 ) 

STRAGE ( I )-0. 

STRCHG ( I ) -STRCHG ( I ) -WTRDEF ( I ) 

IP 

(STRCHG( I-l) .GT. 0. . AND. RMNDR( I ) . GT. 0. ) GOTO 
118 

STRCHG( I)— STRAGE (I-l ) 

WTRDEF ( I )»RMNDR( I )-STRCHG( I ) 

118 CONTINUE 
C 

C DO LOOP 20 CALCULATES THE YEARLY SUMS 

OF THE INDICATED 
C QUANTITY 

C 

SUMTF-0. 

SUMPIN-0. 

SUMTC-0. 

SUMPCM-0. 

SUMAPE-0. 

SUMDEF-0. 

SUMSPL-0. 

SUMACT-0. 

DO 20 1=1,12 
SUMTF=SUMTF+TDEGP ( I ) 
SUMPIN=SUMPIN+PPTIN ( I ) 
SUMTC«SUMTC-i-TDEGC ( I ) 
SUMPCM-SUMPCM+PPTCM ( I ) 
SUMAPE=SUMAPE+APETCM( I ) 
SUMDEP=SUMDEF+WTRDEF ( I ) 
SUMSPL=SUMSPL4-WTRSPL( I) 

IF 

(PPTCM(I).GT.APETCM(I) .AND.STRCHG(I) .GE.O. ) 

1 AGTEVP(I)»APETCM(I) 

IF 

(APETCM(I).GT.PPTCM(I) .AND.STRCHG(I) .LE.O. ) 

1 ACTEVP(I)»PPTCM(I)+ABS(STRCHG(I) ) 


SUMACT-SUMACT+ACTEVP ( I ) 

20 CONTINUE 
ADEGF«SUMTF/12. 

ADEGC-SUMTC/12. 
DEXARD»(100.*SUMDEF)/SUMAPE 
DEXHUM- ( 1 00 . *SUMSPLj/SUMAPE 

AMSTDX- ((110. *SUMSPL ) + 60 . *SUMDEF ) /SUMAPE 
RUNOFF-SUMPCM-SUMAGT 

C THE ABOVE THREE STATEMENTS CALCULATE THE 
CLIMATIC INDICES 
C OP THORNTHWAITE 
C 

C TRYING TO GET % OF CAPACITY HERE 
DO 456 11-1,12 

456 PCMAX ( I I ) - ( STRAGE (II) /STRMAX ) * 1 00 . 

WRITE(5,21) 

21 FORMAT (7X,' JAN FEB MAR APR 

MAY JUNE JULY 

1 AUG SEPT OCT NOV DEC ' ) 
WRITE(5,22) (TDEGF(KP) ,KP-1 ,12) 

22 PORMATC TDEGF ’,12F6.2) 

WRITE (5,23) (TDEGC(KP) ,KP«1 ,12) 

23 PORMATC TDEGC ',12F6.2) 

WRITE (5,24) (PPTIN(KP) ,KP=1 ,12) 

24 PORMATC PPTIN ',12P6.2) 

WRITE (5,25) (PPTCM(KP) ,KP»1 ,12) 

25 FORMATC PPTCM ',12P6.2) 

WRITE (5,26) (SMALLI(KP) ,KP-1,12) 

26 FORMATC SMALLI ' , 1 2F6. 2 ) 

WRITE (5,27) (UPETCM(KP) ,KP«1 ,12) 

2 7 PORMATC UPETCM* ,12F6. 2) 

WRITE (5,28) (CORPHI(KP) ,KP»1,12) 

28 PORMATC CORPHI’ ,12F6.2) 

WRITE (5,29) (APETCM(KP) ,KP-1 ,12) 

29 PORMATC APETCM' ,12P6. 2) 

WRITE (5,30) (RMNDR(KP) ,KP»1 ,12) 

30 PORMATC RMNDR ',12F6.2) 

WRITE (5,31) (STRCHG(KP) ,KP-1 ,12) 

31 FORMATC STRCHG' ,12F6. 2) 

write (5,32) (STRAGE(KP) ,KP=1 ,12) 

32 FORMATC STRAGE ', 1 2F6. 2 ) 

WRITE(5,567) ( PCMAX( KP) , KP=1 , 1 2) 

567 FORMATC %CAP ',12P6.0) 

WRITE (5,33) (WTEDEP(KP) ,KP»1 ,12) 

33 FORMATC WTRDEF ' , 1 2P6. 2 ) 

WRITE (5,34) (WTRSPL(KP),KP=1,12) 

34 PORMATC WTRSPL' ,12F6. 2) 

WRITE (5,35) (ACTEVP(KP) ,KP=1,12) 

35 FORMATC ACTEVP ' , 1 2P6. 2 ) 

WRITE(5, 36) 

36 FORMATC 0 ', 22X, 'ANNUAL SUMMARY') 
WRITE(5,37) 
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37 P0RMAT(2X, 'TDEGP TDEGC PPTIN 
PPTCM APETCM WTRDEP 

IWTRSPL ACTEVP') 

WRITE(5,38) 

ADEGP , ADEGC , SUMPIN , SUMPCM , SUMAPE , SUMDEP , 
1SUMSPL,SUMACT 

38 

P0RMAT(1XP6. 2,2X,P6.2,2X,P6.2,2X,P6. 2,2X,P6. 2, 
2X,P6.2 

1,2X,P6.2,2X,P6.2) 


WRITE(5, 39)AMSTDX,DEXARDrDEXHUM,RUNOPP 
39 PORMAT( ' 0 ' , ' THE POUR THORNTHWAITE 

CLIMATIC INDICES: MOISTURE' 

1,P7.2,' DRYNESS' ,P7. 2, ' 

HUMIDITY' ,P7. 2, ' RUNOPP ' ,P7. 2 ) 
IYEARS-IYEARS+1 
102 CONTINUE 

CLOSE (UNIT=21) 

GO TO 962 

STOP 

END 




